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TECHNICAL  NOTE  NO.  1832 

SMALL  BENDING  AND  STRETCHING  OF  SANDWICH-TYRE  SHELLS 
By  Eric  Reissner 


SUMMARY' 


A  theory  has  been  developed  for  small  bending  and  stretching  of 
sandwich-type  shells.  This  theory  is  an  extension  of  the  known  theory 
of  homogeneous  thin  elastic  shells.  It  was  found  that  two  effects  are 
important  in  the  present  problem,  which  have  not  been  considered 
previously  in  the  theory  of  curved  shells:  (l)  The  effect  of  transverse 
Bhear  deformation  and  (2)  the  effect  of  transverse  normal  stress 
deformation.  The  first  of  these  two  effects  has  been  known  to  be  of 
importance  in  the  theory  of  plates  and  beams.  The  second  effect  was 
found  to  occur  in  a  manner  which  is  typical  for  shells  and  has  no 
counterpart  in  flat-plate  theory. 

The  general  results  of  this  report  have  been  applied  to  the 
solution  of  problems  concerning  flat  plates,  circular  rings,  circular 
cylindrical  shells,  and  spherical  shells.  In  each  case  numerical 
examples  have  been  given,  illustrating  the  magnitude  of  the  effects 
of  transverse  shear  and  normal  stress  deformation. 

The  results  of  this  investigation  indicate  the  necessity  of 
taking  account  of  transverse  shear  and  normal  stress  in  sandwich-type 
shells,  as  soon  as  there  is  an  order-of-^nagnitude  difference  between 
the  elastic  constants  of  the  core  layer  and  of  the  face  layers  of  the 
composite  shell.  It  was  found  that  the  changes  due  to  transverse 
shear  and  normal  stress  deformation  in  the  core  may  be  so  large  as  to 
be  no  mere  corrections  to  the  results  of  the  theory  without  transverse 
core  flexibility. 

The  actual  magnitude  of  the  changes  is  greatly  dependent  on  the 
geometry  and  loading  condition  of  the  structure  under  consideration 
so  that  no  general  rules  may  be  given  which  indicate  for  which  elastic 
modulus  ratio  the  changes  begin  to  be  significant. 

Solutions  of  problems  in  the  present  theory  may  in  general  be 
obtained  by  mathematical  methods  which  are  similar  to  those  employed 
in  the  theory  of  plates  and  shells  without  the  effect  of  transverse 
shear  and  normal  stress  deformation  included.  The  present  work  does 
not  include  consideration  of  buckling  and  finite  deflection  effects. 
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INTRODUCTION 


In  this  report  an  extension  of  the  classical  theory  of  small 
"bending  stretching  of  thin  elastic  shells  is  considered.  Instead 
of  a  homogeneous  shell,  investigation  is  made  of  a  shell  constructed 
in  three  layers:  A  core  layer  of  thickness  h  with  elastic 
constants  Ec,  Gc,  and  vc  and  two  face  layers  of  thickness  t  with 

elastic  constants  Ef,  Of,  and  Vf.  In  the  developments  certain 

restrictive  assumptions  are  made  which  somewhat  limit  the  general 
applicability  of  the  results.  In  so  doing  formulas  are  obtained  which 
are  as  compact  as  possible  while  still  describing  the  essential 
characteristics  of  the  sandwich— type  shell. 

The  thickness  ratio  t/h  is  assumed  small  compared  with  unity; 
at  the  same  time  the  ratio  Eft/Ech  is  assumed  large  compared  with 
unity .  This  latter  assumption  means  that  the  face  material  is  so 
much  stiffer  than  the  core  material  that  the  contribution  of  the  core 
layer  to  stress  couples  and  tangential  stress  resultants  of  the  composite 
shell  is  negligible.  It  is  known  that  for  flat  plates  these  assumptions 
necessitate  the  taking  into  account  of  the  effect  of  transverse  shear 
deformation.  (See,  for  instance,  reference  1.)  Th9  same  would  be 
expected  to  be  true  for  curved  shells,  and  the  present  report,  therefore, 
gives  a  system  of  equations  in  which  this  effect  is  incorporated. 

A  further  effect  which,  it  appears,  has  not  been  considered 
previously  in  the  analysis  of  small  deflections  of  sandwich  structures 
is  the  effect  of  transverse  normal  stress  deformation.  In  the  present 
report  it  is  shown  that  this  effect  arises  in  a  manner  which  is  typical 
for  shells  and  has  no  counterpart  in  plate  theory.  It  may  be  likened, 
roughly,  to  what  happens  in  the  bending  of  curved  tubes. 

The  process  by  which  the  general  results  of  this  report  are 
obtained  iB  as  follows:  First,  each  of  the  face  layers  of  thickness  t 
is  assumed  to  behave  like  a  thin  shell  without  bending  stiffness.  The 
loads  applied  to  these  face  shells,  henceforth  called  face  membranes, 
are- of  two  kinds:  (l)  External  loads  and  (2)  loads  caused  by  the  stresses 
in  the  core  layer.  Next,  the  core  layer  of  thickness  h  is  assumed 
to  behave  like  a  three-dimensional  elastic  continuum  in  which  those 
stresses  which  are  parallel  to  the  faces  are  negligible  compared  with 
the  transverse  shear  and  normal  stresses .  On  the  basis  of  these  two 
assumptions  three  steps  are  carried  out.  First,  the  equilibrium 
equations  of  the  core  layer  and  of  the  face  layers  are  obtained.  Then 
an  appropriate  expression  for  the  strain  energy  of  the  composite 
structure  is  derived.  Finally,  Castigliano*s  theorem  of  minimum 
complementary  energy  is  used  to  obtain  the  relations  which  connect 
stress  resultants  and  couples  of  the  composite  shell  with  the  quantities 
yliich-  describe  the  state  of  deformation  of  the  composite  shell. 
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The  system  of  equations  which  is  obtained  in  the  foregoing  manner  is 
specialized  for  the  following  cases: 

(1)  Flat  plate 

(2)  Circular  ring 

(3)  Circular  cylindrical  shell 

(4)  Spherical  shell  with  axisymmetrical  deformation 

In  each  case  a  number  of  problems  are  solved  explicitly  and  the 
appreciable  effect  of  transverse  shear  and/or  normal  stress  deformation 
is  illustrated  numerically. 

This  work  was  conducted  at  the  Massachusetts  Institute  of 
Technology  under  the  sponsorship  and  with  the  financial  assistance 
of  the  National  Advisory  Committee  for  Aeronautics. 


SYMBOLS 


h 

t 

*1'  h 
t 

al*  a2 
^1*  Eg 
Nrmru 

wnm2 
I'm*  Pn2 


core— layer  thickness 

face-layer  thickness 

curvilinear  coordinates  on  middle  surface  of 
composite  shell 

distance  coordinate  measured  along  normal  to 
middle  surface  of  shell 

coefficients  of  linear  element  on  middle  surface 
of  shell 

principal  radii  of  curvature  of  middle  surface  of 
shell 

direct  stress  resultants  in  upper  face  membrane; 
n  =  1,2;  m  =  1,2 

direct  stress  resultants  in  lover  face  membrane 

tangential,  components  of  external  load  intensity 
on  upper  and  lover  membranes 

normal  components  of  external  load  intensity  on 
upper  and  lover  membranes 


k 
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Tl^  T2i 
% 

Tn£u5  Tn(U 

Tn£m 

§1, 

V 

Mrrn 

Pn 

<L 

B 

It 

Ef,  Gf,  V 
Ec>  Gc 

Ml,  U2 

W 

Pi,  P2 
e£m 

C*  =  2tEf 


components  of  transverse  shear  stress  in  core  layer 

component  of  transverse  normal  stress  in  core  layer 

values  of  transverse  shear  stresses  for  £  =  i-  h/2; 
n  =  1,2 

values  of  transverse  normal  stresses  for  £  =  —  h/2 

values  of  transverse  shear  stresses  at  middle  surface 
of  shell 

transverse  shear  stress  resultants 

direct  stress  resultants  parallel  to  middle  surface 
for  composite  shellj  n  =  1,2;  m  =  1,2 

stress  couples  for  composite  shell;  n  =  1,2;  m  =  1,2 

tangential  components  of  external  load  intensity  for 
composite  shell;  n  =  1,2 

normal  component  of  external  load  intensity  for 
composite  shell 

external  load  intensity  term  defined  hy  equation  (22) 

strain  energy 

elastic  moduli  of  isotropic  face-layer  material; 

V  =  Vf 

elastic  moduli  in  transverse  direction  of  core-layer 
material 

effective  tangential  components  of  displacement  of 
elements  of  composite  shell 

effective  normal  component  of  displacement  of  elements 
of  composite  shell 

effective  components  of  change  of  Blope  of  normal 
to  middle  surface  of  composite  shell 

component  of  strain  =  cr^m/Ec^ 


V 


V 
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D*  =  (l/2)t(h  +  t )2Ef 


C  =  C*/(l  -  v2) 
D 

z,  7 

r>  0- 
a 

x,  9 

Xl>  X2>  X12 
4  =  Tt/l 
l 

™2 


■bending  stiffness  factor  (D  =  D*/(l  —  v2)) 

Cartesian  coordinates  in  plane  of  flat  plate 

polar  coordinates  in  plane  of  flat  plate 

radius  of  circular  ring,  cylindrical  shell, 
spherical  shell 

surface  coordinates  on  cylindrical  shell 
parameters  defined  by  equation  (63) 

half  wave  length  of  sinusoidal  load  distribution 
quantities  defined  by  equation  (197) 


k 

0,  9 

(a 

K 


complex  quantity  defined  by  equation  (200) 
surface  coordinates  on  spherical  shell 
quantity  defined  by  equation  (74) 
parameter  defined  by  equation  (190) 


I  —  GENERAL  ^THEORY 
Statics  of  Sandwich-Type  Shell 


In  order  to  derive  a  complete  system  of  equations  for  the  shall 
composed  of  face  layers  and  core  layers  it  is  necessary  first  to 
consider  separately  the  statics  of  the  face  layers  and  that  of  the 
core  layer  of  the  shell.  Combination  of  the  results  obtained  for 
these  two  components  of  the  composite  structure  must  and  will  lead 
to  those  differential  equations  of  equilibrium  which  hold  for  elements 
of  a  shell,  whether  this  Bhell  is  of  homogeneous  or  nonhomogeneous 
construction.  In  addition,  however,  relations  are  obtained  which 
are  characteristic  of  the  sandwich-type  shell. 

Coordinate  system  on  shell. —  A  curvilinear  coordinate  system 
(tq,  1 2*  t)  is  introduced  as  follows:  Let  ^  and  |2  be  coordinates 
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on  the  middle  surface  of  the  composite  shell  and  let  £  be  the  distance 
of  a  point  of  the  shell  from  its  middle  surface,  measured  along  the 
normal  to  the  middle  surface.  In  order  that  this  system  of  coordinates 
be  an  orthogonal  system,  choose  the  |q,  t2  curves  as  lines  of  curvature 

on  the  middle  surface  (in  the  case  of  shells  of  revolution  the  lines  of 
curvature  are  identical  with  the  meridians  and  parallels  on  the  middle 
surface) . 

The  linear  element  in  the  foregoing  system  of  coordinates  is  of 
the  form 


where  and  09  are  the  coefficients  of  the  linear  element  on  the 

middle  surface  and  Rq  and  R2  are  the  principal  radii  of  curvature 
of  the  middle  surface  (see  fig.  l) .  Formulas  for  the  calculation  of 
the  quantities  <%  and  Rn  are  contained  in  texts  on  differential 
geometry.  They  are  collected,  together  with  other  results,  in 
reference  2,  which  deals  with  the  theory  of  homogeneous  thin  shells. 

Statics  of  face  layers.—  The  face  layers  are  treated  as  thin  shells 
of  thickness  t  and  it  is  assumed  that  the  bending  stiffness  of  these 
thin  shells  about  their  own  middle  surface  may  be  neglected.1  Because  of 
this  neglect  from  now  on  they  will  be  designated  as  face  membranes. 

The  middle  surfaces  of  the  face  membranes  evidently  are  given  with 
reference  to  the  three-dimensional  system  of  curvilinear  coordinates 

by  £  =  ^(h  +  t)  and  £  =  —  ?j(h  +  t).  From  equation  (l)  it  follows 

that  the  linear  element  on  the  mdHdle  surfaces  of  the  face  membranes 
is  given  by 


(2) 


The  components  of  external  load  intensity  on  the  upper  and  lower 
membranes  are  designated  by  plu,  p^,  and  and  by  p^,  P22*  011(1 

respectively  (fig.  2) .  The  core-layer  stresses  which  act  on  the  upper 
and  lower  membranes  are  given  as  Tq£u,  t2^u,  and  cx£u  and  by  TgZl* 


1This,  of  course,  means  that  no  local  buckling  phenomena  are 
considered  in  the  present  work. 
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and  (T^  j }  respectively.  Finally,  the  direct  stress  resultants  in  the 
upper  and  lower  face  membranes  are  designated  "by  Nllu,  N12u,  N21u, 
and  W22u  and  "by  N^,  N122i  %L2j  6111,1  n222j  respectively  (fig.  2). 

There  are  then  three  equations  of  force  equilibrium  for  the 
elements  of  each  of  the  two  membranes.  Writing 


onu  -  “11(1  + 
On!  -  “n^  “ 


the  equations  for  the  upper-face  membrane  are  the  following  : 

^SuHUe.  +  +  Nl2u  |S!.  _  „22u  +  -  t1?u)  =  0  (4) 


Sg^1<1^  +  a°^’jllggu  +  N21u  ISu  -  rjr11  +  «iua2u(p2u  -  t2^u)  ■  0  (5) 


”llu  .  ”22ii 


""  4u  +  a£u  “  0 


The  corresponding  equations  for  the  lower— face  membrane  are 


^22%1I  +  g2L^2U  +  wn07  -  H. 


8a- 


222  +  “ll^Z  (p12  +  t1?2)  -  0  (T) 


+  ~ +  %L2  ~  N11Z  +  a12a2l(P22  +  T2Cl)  =  0 


2 These  are  obtained  from  the  corresponding  equations  of  reference  2 
with  On  changed  to  ct^  and  with  stress  couples  and  transverse  shear 
stress  resultants  omitted.  To  make  up  for  this  omission,  the  loads  on 
the  two  membranes  are  assumed  to  act  at  their  middle  surfaces;  this  means 
terms  of  the  order  t/R  are  neglected  (but  not  terms  of  order  h/R) . 


As  'bending  moments  and  transverse  shears  are  assumed  not  to  be 
acting  in  the  individual  membranes  the  moment  equilibrium  equations 
become  the  symmetry  relations 


W12u  =  %lu 
n12 l  =  N21 l 


(10) 


Before  analyzing  the  state  of  stress  in  the  core  layer  it  is 
convenient  to  see  what  relations  follow  from  equations  (4)  to  (9)  for 
the  composite  shell. 

Statics  of  composite  shell.—  It  may  he  seen  that,  in  view  of  the 
fact  that  all  face-parallel  stresses  in  the  core  layer  are  neglected, 
the  following  expressions  for  the  face— parallel  stress  resultants  and  g 

couples  of  the  composite  shell  are  obtained: 
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In  the  same  way  the  following  expressions  are  obtained  for  components 
of  external  force  and  moment  intensity: 


Further ,  a  load  term  of  the  following  form  will  he  encountered: 


h  +  t 


X1- 


(22) 


which  hears  a  relation  to  equation  (20)  similar  to  that  which  equation  (21) 
hears  to  equation  (19) .  This  last  term  would  represent,  for  a  homogeneous 
shelly  the  average  transverse  normal  stress  at  any  station  of  the  shell, 
assuming  that  the  loads  qu  and  q^  alone  are  responsible  for  this  stress. 
For  a  homogeneous  isotropic  shell  this  term  is  of  no  importance.  For  a 
sandwich— type  shell,  as  will  he  seen,  it  may  sometimes  he  of  importance. 


In  order  to  obtain  force  and  moment  equilibrium  equations  for  the 
composite  shell  the  face— membrane  equilibrium  equations  (4)  to  (9)  are 
combined  suitably.  Adding  equations  (4)  and  (7),  and  (5)  and  (8), 
respectively,  the  two  equilibrium  equations  for  the  force  components 
parallel  to  the  middle  surface  of  the  shell  are  obtained.  In  order  to 


reduce  them  to  known  form  (see  reference  2)  the  following  relations  are 

used  between  the  core-layer-surface  shear  stresses  t  e  t  » 

nsu  ntZ 

and  the  transverse  shear  stress  resultants  0-i  and  Cu. 
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Equations  (23)  and  (24)  will  subsequently  be  shown  to  be  in  agreement 
with  the  usual  definition  for  the  transverse  .shear  stress  resultants  by 
consideration  of  the  stress  distribution  of  the  core  layer. 

With  equations  (23)  and  (24) ,  there  are  obtained  by  combination  of 
equations  (4)  and  (7),  and  (5)  and  (8)  -  carrying  out  addition  as  well 
as  subtraction  -  the  following  four  equations: 


(25)* 


vsfr1  +  +  IS _  ^  sir  +  -  %)  - 0 


^0-2^12  8ai^2 

81,  +  Mo 


51?p2  +  %1  |S  -  Ma  |a  +  ajc^  -  %)  -  0 


Two  further  equations  are  obtained  by  adding  and  subtracting, 
respectively,  equations  (6)  and  (9).  Adding  equations  (6)  and  (9)  and 
taking  account  of  equations  (ll),  (l4),  and  (20),  there  follows: 
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In  order  that  this  reduces  to  the  correct  equation  of  transverse  force 
equilibrium  as  given  in  reference  2, 


a^Og 


j1 ♦ wX1 + wh*  -  i1  -  wX1  -  wk 


Sa2*5l  SaiQ2 
dlq  dl2 


(30) 


Equation  (30),  Just  as  equations  (23)  and  (2k),  can  again  be  verified 
independently  "by  consideration  of  the  state  of  stress  in  the  core  layer. 
On  the  basis  of  equation  (30),  equation  (29)  is  written  in  the  form  . 


dagQl 

Bli 


+  a-jOgq. 


0 


(3D* 


The  last  equation,  use  of  which  is  required  for  the  sandwich-type 
shell  and  which  has  not  previously  been  given,  is  obtained  by 
subtracting  equation  (9)  from  equation  (6).  Taking  account  of 
equations  (15),  (l8),  and  (22),  there  results 


2ala2 
h  +  t 


^22 \  _  2ana0s  +  a-,a0  / 

1  +  *L±J*\( 

1  + 

VE1 

e2  /  12  1  2 |_\ 

®2/' 

2El  ) 

-  wX1  ~  w)^] = ° 


os) 


Provisionally ^  there  is  written 


(*  +  +  ^XEu  +  i1  -  '  2"£m  ( 33) 


V 
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and  it  will  subsequently  be  shown  that  represents  the  value 

of  cr£  at  the  middle  surface  of  the  shell.  Combining  equations  (33) 
and  (32)  yields 


+ 


1  (* 11 
h  +  t\E^ 


s  =  0 


(3b)* 


Equation  (34)  has  no  relation  to  the  sixth  equation  of  equilibrium 
for  an  element  of  the  shell  which  expresses  the  condition  of  moment 
equilibrium  about  the  normal  to  the  middle  surface.  That  equation  which, 
as  is  known,  is  an  identity  when  resultants  and  couples  are  expressed 
in  terms  of  stresses  does  not  occur  in  the  present  derivations,  or  rather 
it  is  contained  in  equations  (12),  (13),  (l6),  and  (17),  which  give 
explicitly  the  slight  differences  between  Ni  p  and  Npi  ,  and  M-]  p 

and  M21 . 

Stress  distribution  in  core  layer.—  In  order  to  verify  independently 
equations  (23),  (24),  and  (30),  as  well  as  for  the  subsequent  derivation 
of  appropriate  stress-strain  relations,  it  is  necessary  to  determine  the 
distribution  of  stress  in  the  core  layer. 

Assuming  that  the  components  of  stress  ffp,  and  in  the 

core  which  would  contribute  to  stress  resultants  and  couples  of  the 
composite  shell  are  of  negligible  importance, 3  these  components  of 
stress  may  be  set  equal  to  zero  and  only  the  components  of  transverse 
shear  stress  and  transverse  normal  stress  t^,  t2£  ,  and  erf;  may  be 

retained.  The  differential  equations  of  equilibrium  for  these  three 
remaining  components  of  stress  in  the  system  of  curvilinear  coordinates 
defined  by  equation  (l)  are  obtained,  from  the  general  form  of  these 
differential  equations  in  reference  3,  in  the  following  form: 


3lt  is  for  this  purpose  that  the  order— of-magnitude  relation 
hEc/tEf  «  1  is  assumed. 
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The  values  of  the  three  stress  components  at  the  middle 
surface  ( t  =  0)  are  designated  by  the  subscript  m.  Integration  of 
equations  (35)  to  (37)  then  gives 


ru 


(1  +  i/sxf(i  4  tfe) 


(38) 


The  transverse  shear  stress  resultants  Qx  and  Q2  are  obtained 
from  equations  (38)  and  (39)  in  the  form1' 

H(h+t)/2 

%'L  (hta-rrfa) 


-  + 

"  ,  /h  +  tf 

V2IW 

The  integration  must  be  extended  over  the  thickness  of  the  core 

a}8°  over  the  thickness  of  the  face  layers,  in  accordance 
with  the  prior  assumption  that  the  stresses  Tn£Uj  Tn£z#  a£u, 

and  ofr  may  be  taken  to  act  at  the  middle  surfaces  of  the  respective 
face  membranes . 
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Wow  aB  intended,  the  proof  is  carried  out  of  equations  (23) ,  (24), 
and  (30),  which  were  used  to  obtain  the  differential  equations  for  the 
composite  shell . 

To  verify  equation  (23),  from  equations  (38)  and  (39)  for  the 
left-hand  side  of  equation  (23),  the  following  equation  is  obtained: 


TnU 


1  + 


h  +  t 
2*n 


1  - 


Tn£m 
h  +  t 


2*n 


h  +  t 


Tn£m 


-m 


and  this,  in  conjunction  with  equation  (4l),  verifies  equation  (23). 

To  verify  equation  (24)  in  the  same  manner,  from  equations  (38) 
and  (39)  for  the  left  side  of  equation  (24),  the  following  equation 
is  obtained: 


and  this,  in  conjunction  with  equation  (4l),  verifies  equation  (24). 

To  verify  equation  (30),  equation  (40)  is  used  to  write  for  the 
left  side  of  equation  (30) 


/ «2TlSm  V  8  fa 1t2U  \  5  /a2TlU  \  + 

siiL  +  h+_t ) +  ai2(i  +  2u_±)  sii(i_1l±1 


2R- 


2 *2/ 


^  /  ga2Tl^m 


5  /  galT2U 


hjf  t 

2R, 


2R- 


2  |  +  81 


2!  1  _  fh  +  t 


2R 


and  this,  in  conjunction  with  equation  (4l),  verifies  equation  (30) 
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The  Bection  on  "the  stress  distribution  in  the  core  layer  is 
concluded  hy  listing  the  form  which  equations  (38)  to  (40)  for  the 
stresses  in  the  core  layers  assume  for  "thin"  shells,  that  is  for 
shells  for  which  h/R  «  1.  From  equations  (38)  and  (39)  in* 
conjunction  with  equation  (41),  it  follows  that 


Q-, 

Tl£  “  h  +  t 


To/-  =  Qg 

2£  h  +  t 


m* 


From  equation  (40),  in  conjunction  with  equation  (4l),  it  follows  that 


“l^m  ~ 


h  + 


+ 


(43a) 


It  is  necessary  to  note  for  some  of  the  following  considerations 
that,  in  view  of  equation  (31),  instead  of  equation  (43a)  there  may 
be  written 


at  =  +  If  - *)  (43b)* 


It  is  seen  that  in  this  approximation  the  transverse  shear  stresses  are 
uniform  across  the.  thickness  of  the  core  layer,  while  the  transverse 
normal  stress  is  composed  of  two  terms,  one  uniform  across  the  thickness 
and  the  other  varying  linearly  across  the  thickness. 

No  further  calculations  are  needed  with  reference  to  the  state 
of  stress  in  the  composite  shell.  The  next  step  is  to  complete  the 
system  of  differential  equations  for  stress  resultants  and  couples  by 
deriving  an  appropriate  system  of  stress-strain  relations. 


Strain  Energy  of  Sandwich-Type  Shell 

In  calculating  the  strain  energy  of  face  membranes  and  core  layer 
it  is  assumed  that  both  are  isotropic  and  elastic,  with  elastic 
constants  Ef,  Vf  =  V,  Gf  =  Ef/2(l  +  V )  and  Ec,  Vc,  Gc  =  Ec/2(l  +  Vc). 
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Poisson's  ratio  for  the  face  membranes  is  written  without  a  subscript, 
because,  in  view  of  the  assumed  stress  distribution,  there  is  no  explicit 
occurrence  of  Poisson's  ratio  vc  for  the  core  layer. 

The  strain  energy  for  the  composite  shell  is  the  sum  of  the  strain 
energies  for  the  face  membranes  and  for  the  core  layer 


it  -  Tif 


For  the  purpose  of  obtaining  stress-strain  relations,  both  itf  and  jtc 
are  expressed  in  terms  of  stresses  rather  than  in  terms  of  strains. 

Strain  energy  of  face  layers.-  Considering  that  the  element  of 
area  on  the  middle  surfaces  of  the  membranes  is  of  the 
form  Wl  ±  ^)(l  ±  «1  and  that  the  stresses 

in  the  membranes  are  the  stress  resultants  divided  ^  the  membrane 
thickness  t,  there  is,  from  well-known  principles,  the  following 
relation: 


■  tJJ  + 1,221,2  -  2vwn“i,2a‘ + 2(1 + v)Hi2u2] 
x  (a + WX1  +  WX1"2 dSi  162 

+  iff *  "22T2  -  2V%11N222  +  2<1  +  V)I»12I2] 


*  (l  -  J^)(1  -  W)ai“2  161  312 


Equation  (15)  is  transformed  into  an  expression  containing  stress 
resultants  and  couples  of  the  composite  shell  by  means  of  equations  (11) 
to  (18)  which  lead  to  the  relations 
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/ 

« 


1 


2(l  +  ¥)%.  ’  *11  +  rh  «L1 
2(1  -  -  »11  -  rfr  %1 


(46) 


with  corresponding  formulas  for  N-]  p  and  ^22  ^ 

In  this  transformation  the  cases  are  limited  to  those  for 
which  h/R  « 1.  Then,  with  the  two  constants  C*  and  D*  defined  by 


C*  =  2tEf 


D*  =  1  t(h  +  tJ^Ef 


the  following  expression  for  itf  is  obtained: 


*f  =  2 


in 

Ji 


r- 

-  Jj[nu2  +  n22?  -  2vKui(22  +  2(1  +  vlHj/J 


+  d?K2  +  1^2^  —  2vMjj_Mgg'  "t  2(l  +  v)M22^||'  ®2®2  ^^1  ^-^2  (48)* 


It  may  be  remarked  that  equation  (48)  could  have  been  given  directly,  by 
analogy  with  known  results  for  the  isotropic  homogeneous  shell. 

Strain  energy  of  core  layer.—  With  the  stresses  o±,  cr2,  Ti2 

assumed  to  vanish,  there  results  for  the  strain  energy  of  the  core  layer 


5lTote  that  equations  (46)  and  corresponding  equations  can  be  used 
to  calculate  the  stresses  in  the  two  different  face  membranes,  once 
stress  resultants  and  couples  in  the  composite  shell  are  known. 


4 
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jr  =2 
c  2 


+  T2l 


('*«(■ 


d^a. 


Again  the  terms  £/R  compared  with  unity  are  neglected  and,  consistent 
with  this  neglect,  the  values  of  the  stresses  Tn£  and  cr£  are  taken 

from  equations  (42)  and  (43) • 

The  value  of  may  he  chosen  from  either  equation  (43a)  or 

equation  (43b).  The  form  of  the  results  depends  somewhat  on  which  of 
the  two  equations  is  chosen,  in  the  sense  that  the  meaning  of  the 
deformation  quantities  which  are  to  he  determined  depends  on  which 
of  the  two  equations  is  taken.  This  question  is  decided  in  the 
following  manner:  As  all  resultants  and  couples  enter  the  expression 
for  the  strain  energy  only  as  themselves  and  not  in  differentiated 
form,  except  when  equation  (43a)  is  used,  the  selection  of  equation  (43b) 
for  is  proposed,  thereby  excluding  derivatives  of  stress  resultants 

and  couples  from  the  expression  for  the  strain  energy  it. 

Introducing  then  equation  (43b)  into  equation  (49)  yields 


«c  =  2 


Qi2  +  Q22  _L 

(h  +  t)2cc  E° 


[°u  -  *  s?  - «)]  ' 


X  df^aqC^  dlq  d^ 


The  integration  with  respect  to  £  is  carried  out  and 
equation  (50)  becomes 


*c  =  2 


Qq2  +  Qe2 
7FTTTgT 


h  +  t 


aU2  + 


1  /«n 

12U 


It  was  to  be  expected  that  the  terms  containing  the  modulus  of 
rigidity  Gc  would  occur  in  the  foregoing  form.  The  contribution  of 
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the  present  report  up  to  this  point,  besides  giving  the  new  equation  (34) 
for  cr£m,  is  thought  to  be  the  determination  of  the  form  in  which  the 

effect  of  transverse  normal  stress  deformability  manifests  itself  in 
the  strain  energy  of  the  sandwich  shell. 


Stress-Strain  Relations  for  Composite  Shell 

In  what  follows  a  system  of  stress— strain  relations  for  the  com¬ 
posite  shell  is  obtained  by  the  use  of  Castigliano’s  theorem  of  minimum 
complementary  energy.  The  manner  in  which  the  theorem  is  used  here 
appears  to  have  been  employed  first  by  E.  Trefftz  (reference  4)  for  the 
purpose  of  avoiding  geometrical  considerations  in  the  derivation  of  the 
stress— strain  relations  for  thin  homogeneous  shells  with  smal 1  defor¬ 
mations,  without  consideration  of  the  effects  of  transverse  shear  and 
normal  stress  deformation. 

Assuming  for  the  present  purpose  that  all  boundary  conditions  for 
the  shell  under  consideration  are  stress  conditions,  the  theorem 
consists  in  the  statement  that  among  all  statically  correct  states  of 
stress  the  actually  occurring  state  of  stress  makes  the  strain  energy 
of  the  system  a  minimum.  In  the  application  of  the  theorem  the  fact 
is  taken  into  account  that  statically  correct  states  of  stress  only 
are  to  be  compared,  by  means  of  the  Lagrangian  multiplier  method. 

Before  minimizing  *  an  integral  is  added  to  it  which  contains  the 
six  equilibrium  equations  (25)  to  (28),  (31),  and  (34),  each  of  the 
six  equations  multiplied  by  a  Lagrangian  multiplier.  It  can  then 
be  shown,  by  using  Castigliano’s  theorem  with  prescribed  boundary 
displacements  instead  of  with  prescribed  boundary  stresses,  that 
each  of  the  six  multipliers  has  the  meaning  of  one  of  the  displacement 
quantities  which  occur  in  the  shell  problem. 6 

With  the  foregoing  understanding  of  the  meaning  of  the  multipliers, 
the  multiplier  of  equation  (25)  is  designated  by  uq;  that  of 

equation  (26),  by  U2;  that  of  equation  (27),  by  pq;  that  of 

equation  (28),  by  02;  that  of  equation  (31),  by  wj  and  finally  that 

of  equation  (34),  by  k.  It  is  known  that  uq,  u2,  and  w  represent 
the  effective  components  of  displacement  in  the  ij,  §2,  and  £  direc¬ 
tions,  respectively.  Further,  it  is  known  that  $q  and  p2  represent 
the  angles  through  which  the  normal  ±0  the  middle  surface  of  the  shell 
turns  toward  the  iq  and  £2  curves,  respectively.  There  is  no 

^For  the  special  case  of  the  flat  plate  this  haB  been  carried  out 
explicitly  in  reference  1.  For  the  case  of  the  homogeneous  shell, 
without  effect  of  transverse  shear  and  normal  stress  deformation,  the 
proof  has  been  given  in  reference  4.  The  proof  for  the  more  general 
case  which  is  here  considered  is  not  included  as  it  does  not  offer  any 
clearer  insight  into  the  problem  and  tends  to  lengthen  the  analytical 
discussion. 
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immediate  simple  geometrical  interpretation  for  k  and,  while  such 
interpretation  in  terms  of  an  average  transverse  normal  strain  might 
he  deduced  herein,  k  is  considered  as  an  auxiliary  variable  presently 
to  he  eliminated. 

Combining  now  equations  (44),  (48),  (51) ,  and  (25)  to  (28),  (31), 
nnrl  (3k)  in  the  manner  indicated,  the  following  variational  equation 
results: 


no 


k  5 


UU 


^p[%l2  +  N222  -  2VW11W22  +  2(1  +  V)W122] 

+  +  ^22  “  2VML1^2  +  2(1  +  V)M122] 


+  %2  +  Q22  +  h  +  t 
+  (h  +  t)Gc  Ec  L 

m 

+  8  1 1  ^u;l 


>0,30,2  d£-j_ 


Uli 


ScteNn  .  SaiN2i  ,  w  8ai  m  80,2 

^  ~  +  Nt  +  N12  ^7 - w22  5— 


8i 


81 


81, 


bh 


/Q 1 

+  aia2bj  +  pi 


da2%2  ^1^22 

8|j  8S2 


802 


+  ^2 11  ^  —  N- 


8a- 


21  8?7  ”  "H  81 2 

*  *)] ♦  ^  + 


8a2 


8ax 


+  ^  -  %)] ♦ ^  sj  -  %i  r 
+  -  <te)]  +  +  "  “l"^1  +  sf2  ~  « 
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The  variations  in  equation  (52)  are  carried  out  and  integration  is 
done  by  parts  to  eliminate  derivatives  of  variations  in  the  double 
integral.  The  line  integrals  along  the  boundary  which  occur  due  to 
this  integration  by  parts  vanish,  because  it  has  been  assumed  that  all 
stresses  are  prescribed  at  the  boundary  and  therefore  their  variations 
vanish  at  the  boundary. 

The  resultant  variational  equation  is 


nn 


dd 


r 


fell  -  VN22  _  _L  _  JJ2_  _  v.  .  h  +  t  All  ,  N22  A 
°  11  [  C*  al  ala2  Si2  El  +  12EcEi\£i  +  B2  7 


.  bit  TN22  ~  VIJ11  1  ^u2  ui  Sa£  _  _w 

22|_  c*  ®2  5£2  al“2  d£q  Eg 


4.  h  +  t  pil  ,  ^22  _  n 

»  i  -n  I  n  *t  -q  y 

E2  / 


12ECE2\E^ 


+  8N12 


u^  5a^  ^2  u2  2 

+  5|^  ~  al  +  ^2 


2(1  +  V)Ui2  1  5uj 

C*  ctg  dl2 


+  SMup^1  w-2^ 


5Px  P2  5a 


V_A^L*| 
^  + 


k 

cx^cc^ 


al  5ix  al“2  51 2  (h  +  t)Ei  “l0^ 

+  6Mpof^g.  ~  VMl1  -  i  _  _Jl_  ^02 
E*  «2  5I2  0x02  ^ 

'  +  «J2!^  V,M 


(h  +  t)R2  ccjctQ 


12  1  dPi 

«2  dl2 


% 

_(h  +  ' 


t)G~  +  R 


,  Pi  Saa  _  _jl_  5p£  ,  P2  <W]  .  cq 
«1“2  5|2  al  51 1  +  al“2  51  jJ  + 


UI 

1 


V-1  ®{m  +  ;rV 

iiiQ  QC-^OC^ 


>  aia2  d d i2  =0 


(53) 
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As  all  nine  variations  in  equation  (53)  are  independent  of  each  other, 
it  follows  that  the  contents  of  all  nine  brackets  in  equation  (53)  rauBt 
vanish  separately.  Thus  the  following  nine  stress-strain  relations  are 
obtained  for  the  sandwich  shell,  indicating  with  an  asterisk  those 
which  appear  in  final  form. 


W-  (h  +  t)(?A  _  ^22  L  _  (h  1  t)C*\ 

C*  V  VS&f)  C*  V  12Ec%E2  ) 


=  jl— i  +  +  jsl  +  (h  +  Mi 

ai  51-l  +  0302  d|2  E1  12EcEi 


(fa  ±  t)cA  _  WV  _ 

i2ece22  y  C*  \ 


(h  +  t)cA 

12EcBjE2  y 


1  ui  ^2  y  (h  +  t)q 

a2  SI2  +  ala2  Siq  +  E2  12EcE2 


2(1  + 
c* 

—  n12  = 

Ol 

02  MgW; 

M].i  —  VMpp 

=  _L 

Ml 

.  _j£_ 

D* 

al 

aii 

ala2  Sl2 

5? 

1 

CVJ 

£ 

=  -L 

832 

+  -h-  Ms 

D* 

“2 

Sl2 

“l0^ 

02  d  fa 2\ 

ai  »iW 


“l0^  (h  +  t)Ex 


_ _ k _ 1 _ _ 

a,l<X2  (h  +  t)B2 


2(1  +  V)  „  =  ai  5  /Pj\  +  o£  5  /P2\ 

D*  12  02  bi2\alJ  al 


m 

(55) 

(56) * 

(57) 

(58) 

(59) * 
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(h  +  t)Gc 


Pi 


+  J=_ 

“1  Six  % 


(60)* 


% 

(h  +  t)Gc 


^2 

E2 


(61)* 


(h  +  t) 


(62) 


It  may  be  verified  that  the  meaning  of  the  quantities  u2,  w,  0^ 

and  02  is  as  ^as  ^aan  indicated  by  comparing  equations  (5*0  to  (6l) 

with  the  corresponding  equations  of  reference  2  for  the  homogeneous 
shell  with  Ec  =  G-c  =  00. 

The  system  of  equations  (5*0  to  (62)  may  be  brought  into  a 
slightly  more  concise  form  as  follows:  Define  the  quantities  X,2, 

and  ^12  "by 


11 


(h 


t)t  Ef 
El  2  Ec 


\2 


1  (h  +  t)t 

2  ~R^  Ec 


Xl2  = 


1  (h  +  t)t  Ef 

2  RiEg  Ec 


(63) 


and  eliminate  k  from  equations  (57)  and  (58)  by  means  of  equation  (62) 
and  the  equilibrium  equation  (34).  Retain  equations  (56)  and  (59)  to  (6l) 
in  the  foregoing  form  and  write  for  equations  (54)  and  (55) 


(l+l  0Nl1  “  (v  ■  3  Xl2)N22  = 


U2 

aja^ 


dai 


(h  +  t)C* 
12ECR! 


4) 


(64) 
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^1  +  i  X^)n22  “  (v  "  3  X12^)N11 


Equations  (57)  and  (58)  become,  if  D*/(h  +  t)Ec  =  i  t(h  +  t)Ef/Ec, 
according  to  equation  (4-7), 


(l  +  Vl)Mli  -  (v  -  Xi2)^2  =  +.  c|%  |f^  +  ^  8  (66>* 


(1  +  -  (V  -  M.2)Mii  =  |f|  +  ^  ^  ^  s  (67)* 

With  these,  last  transformations  there  is  obtained  a  system  of  V 

equations  which  is  formally  equivalent  to  the  corresponding  system  of 
equations  for  the  homogeneous  shell.  The  5  equilibrium  equations  (25) 
to  (28)  and  (31)  and  the  8  stress-strain  relations  (5 6),  (59),  (60),  (6l), 
and  (6k)  to  (67)  are  used  for  the  determination  of  13  quantities:  Five 
stress  resultants  N-q,  N22,  Q-j_,  and  Qq;  three  stress  couples  Mjj_, 

M,2,  and  Mi2>  and  five  displacements  and  changes  of  slope  u^,  u^  w,  P-^, 
and  P2.  The  quantity  cr(;m  which  occurs  in  the  sixth  equilibrium 
equation  (equation  (34) )  may  be  determined  directly,  once  the  shell 
bending  and  stretching  problem  has  been  solved. 

It  is  seen  that  the  effect  of  transverse  shear  deformation  enters 
equations  (60)  and  (6l)  only  and  that,  when  Gc  =  00,  these  equations 
give  the  values  of  the  known  theory  of  homogeneous  shells  without 
transverse  shear  deformation  (references  2,  3,  and  4) . 

The  effect  of  transverse  normal  stress  deformation  enters  equa¬ 
tions  (64)  to  (67)  only.  It  is  seen  that  it  is,  in  part,  responsible 
for  the  occurrence  of  apparent  stiffness  factors  C*/(l  +  X) 
and  D*/(l  +  X).  Thus,  according  to  equation  (63),  the  effect  of 

finite  Ec  is  to  make  the  shell  more  flexible  in  bending  and  stretching  « 

than  it  would  be  with  Ec  =  «>.  This  effect,  however,  is  present  only 
in  curved  structures  and  not  in  plates  and  straight  beams,  as  the 

quantities  X  have  one  or  both  of  the  radii  of  curvature  in  the  * 

denominator.  A  further  effect  of  finite  Ec  is  occurrence  of  the 


n*/jL  .  U1  d<i2  , 

V  I  /V  V  »  T  /-V  /-V  V  k  » 


°2  al“2  dSq  E2 


2-) 

R2) 


(h  +  t)C* 

+  - - —  q 

12EcR2 


(65) 
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external  load  terms  q  and  s  in  the  stress-strain  relations.  Both 
these  effects  represent,  roughly  speaking,  what  happens  to  the  shape 
of  an  element  of  the  composite  shell  if  the  length  of  the  core  filers 
in  transverse  direction  is  changed,  without  any  stretching  or 
compressing  of  the  face-membrane  elements. 

Having  derived  the  general  system  of  equations  for  the  small 
lending  stretching  of  sandwich— type  shells ,  it  remains  to  apply 
these  equations  to  specific  problems  which  may  le  of  interest  and  to 
determine  the  quantitative  effect  of  the  terms  which  are  characteristic 
of  the  sandwich-type  Bhell.  Some  of  this  work  is  done  in  part  II  of 
the  present  report,  which  follows. 

It  may  le  stated  once  more  that  for  these  specific  applications 
the  five  equilibrium  equations  (25)  to  (28)  and  (31)  and  the  eight 
stress-strain  relations  (56),  (59),  (60),  (6l),  and  (64)  to  (67)  are 
used. 


II  -  APPLICATIONS  OF  GENERAL  THEORY 
Flat  Plates 


The  problem  of  the  flat  plate  is  considered  first  in  order  to  show 
that  the  results  of  reference  1  are  contained  in  the  present  results 
and  in  order  to  solve  some  problems  in  the  theory  of  plates  which  have 
not  been  solved  in  reference  1. 

Rectangular  plates.-  Using  notation  which  is  customary  in  plate 
theory  there  is  set 


ll  =  X 

12  =  Y 

<X]_  =  0C2  -  1 

E1  =  E2 

3 

li 

H 

2 

Ug  =  V 

Pl  =  Px 

P2  =  Py 

Nil  =  % 

>> 

II 

CVJ 

eS"1 

II 

OJ 

%  =  9x 

K? 

II 

<© 

11 

rH 

& 

m12  =  **ry 

* 

II 

OJ 

Pi  =  Px 

P2  =  Py 

ml  =  “x 

m2  =  my 

(68) 
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The  equilibrium  equations  (25)  to  (28)  and  (31)  become 


^Nyy  BNy 


|Si  +  !i  +  q.o 

OX  OJ 


dMg  ,  dMxy 
cfcc  +  c)y 


-  Qx  +  %  = 


+  _  q  +  n,  = 

bx  by  ^ 


The  stress— strain  relations  (5 6),  (59) ,  ( 60) ,  (6l),  and  (64) 
to  (67)  become 


N  -  VNV  =  C*  ^ 
*  J  bx 


W7  -  VI*I  =  °*  | 


2(1  +  VlN^  =  C*(|  ♦  &) 


Sr  -  <h  +  t)<5c(fi  +  ||) 


Qy  =  (h  +  t)GC(Py  +  V- 
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Mjj-  -  VM^  =  D* 


My-  —  VMj  =  D* 


8x 

BiPj 

by 


2(1  +  V 


dy  +  bx 


(73) 


As  in  the  small-deflection  theory  of  homogeneous  plates,  the 
equations  for  stretching  (equations  (69)  and  (71))  are  independent  of 
the  remaining  equations  for  transverse  bending.  Equations  (69)  and  (71) 
for  the  stretching  are  not  affected  by  the  elastic  properties  of  the 
core  layers . 

Equations  (70)*  (72),  and  (73)  have  been  treated  in  reference  1  by 
means  of  a  stress  function  t,  which,  together  with  the  deflection  w, 
was  taken  as  one  of  two  basic  variables.  In  what  follows  an  alternate 
treatment  is  given,  in  which  the  problem  is  reduced  to  three  simultaneous 
equations  for  the  quantities  Px*  Py*  and  w.  On  the  "basis  of  these  three 
simultaneous  equations  a  problem  non  considered  in  reference  1  is  treated, 
namely,  the  bending  of  a  rectangular  plate  which  is  simply  supported  on 
all  four  edges. ? 

To  reduce  equations  (70),  (72),  and  (73)  to  three  simultaneous 
equations  for  $x,  and  w,  first  a  quantity  co  is  defined  by 


dy 


(7*0 


Introducing  equation  (72)  into  the  first  of  equations  ( 70) ,  in  view  of 
equation  (7^),  there  is  obtained 


tt>  +  V2W  =  -  q/(h  +  t)Gc  (75) 

Next,  Qx,  Mjj-,  and  are  taken  from  equations  (72)  and  (73)  and  the 

result  is  substituted  in  the  second  of  equations  (70).  This  gives,  after 
slight  transformations. 


7 

This  same  problem  has  also  been  solved  by  L.  H.  Donnell  by  a 
method  which  differs  from  the  one  employed  here.  (See  reference  5 
where  the  case  of  the  homogeneous  plate  is  considered.) 
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5£1_  V2|3x  _  2(h  +  t)Gc3z  +  -  2(h  +  t)Gcw]  +  %  =  0  (7 6) 


In  an  analogous  manner  the  following  further  equation  is  obtained: 


v2 3y  -  2(h  +  t)Gcpy  +  ^\f^v  ~  2(h  +  t)GcW]  +  “y  =  °  (77) 


In  order  to  solve  equations  (75)  to  (77)  two  equations  are  next  obtained 
involving  w  and  o>  only.  Differentiating  equation  (76)  with  respect 
to  x  and  equation  (77)  with  respect  to  y  and  adding  the  two  resultant 
equations,  in  view  of  equation  (7^),  gives 


2(h  +  t)Gc(co  +  V2w)+|^  +  ^  = 


and,  making  use  of  equation  (75), 


The  following  procedure  may  now  be  carried  out:  (a)  Solve 
equation  (78)  for  o>,  (b)  with  this  value  of  solve  equation  (75) 
for  w,  (c)  substitute  and  w  in  equations  (76)  and  (77)  and  solve 
for  3X  and  p y,  and  (d)  elimipate  extraneous  terms  in  3X  and  Py 

by  considering  equation  (lb) . 

Before  deriving  the  solution  of  a  problem  along  these  lines, 
the  explicit  differential  equation  for  w  which  follows  by  combining 
equations  (75)  and  (78)  may  be  given 


Note  that  the  effect  of  transverse  shear  occurs  on  the  right  side  of 
the  equation  only.  In  order  to  compare  the  magnitude  of  the  q  terms 
on  .the  right  of  equation  (79),  assume  that  relevant  changes  of  q 
occur  over  distances  of  order  l  (where  l  may  or  may  not  be  a 
representative  diameter  of  the  plate).  Then,  as  order-of-magnitude 
relations,  there  results 
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Prom  equation  (80),  it  follows  that  transverse  shear  ceases  to  he  a 
secondary  effect  as  soon  as  Z  is  of  order  \Jht  \JSf /G-c  or  of  smaller 
order . 

Bending  of  rectangular  plate  with  simply  supported  edges.-  The  edges 
of  the  plate  are  assumed  to  he  at  x  =  0,a  and  y  =  0,h  and  along  these 
edges  moments  and  deflections  are  assumed  to  vanish.  Further, 


mX  =  Ely  =  0 


-11 

m=l  n=l 


tmn  Bin  \lx  sin  Mn7 


(81) 


vhere 


=  m*/a 
Mn  =  njt/b 


(82) 


From  equation  (78),  it  follows  that 


_ -km _ 

,  2  2 

^m  +  Mn 


sin  sin  uny  +  qh 


where  is  a  harmonic  function.  Putting  equation  (83)  into 

equation  (75), 


(83) 
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v^-II 


^mn 


l^D _  +  _ L 


V2  +  Mn2  (h  +  t)Gc 


sin  sin  ^7  - 


which  is  integrated  to 


-ill 


3-mn 


(\i2  +  Mn2) 


l  D(Xm2  +  Mn2) 

(h  +  t)G0 


sin  sin  j^y  +  wh  (84) 


where  w^  is  the  general  solution  of  V^wfr  =  -0*1.  It  is  to  he  expected 
and  may  he  shown  explicitly  that  for  the  plate  which  is  simply  supported 
all  around  w^  =  =  ^  anci«?  aB  "^e  Wavier  solution  for  the  plate 

without  transverse  shear  deformation,  the  particular  integral  is  the 
complete  solution  of  the  problem. 

Equation  (84)  may  he  rewritten  in  the  more  explicit  form 


w 


Qmn" 

Ef  (h  +  t)t 

1  jjn2  +  n2(a2/b2)j| 

L1  2  .(!• 

-  V2)Gr  a2 

m2  +  n2  (a2/h2)] 

2 

X  sin  Sat  x  sin  ^  y 
a  b 


(85) 


When  Gc  =  oo,  equation  (85)  reduces  to  NavierTs  solution  w^.  Equation  (85) 
is  more  readily  interpreted  hy  means  of  the  ratio  w/wN  of  deflection 

with  and  without  transverse  shear  deformation.  On  the  basis  of 

Q 

equation  (85),  there  may  he  obtained  the  following  equation  (86), 


^Setting  jt^/2(l  —  V2)  =  5.4  and  (Ef/Gc)t(h  +  t)/a2  =  fi, 
equation  (86)  takes  on  a  form  which  contains  as  a  special  case  the 
result  of  equation  (l8)  of  reference  5.  , 
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(h  +  t)gt 
WN 


_  ^  _  (h  +  t)tEjj  V^y-]^ 

"2(l-v2)Gc  ^ 


rt2  Ef 


h  +  t)t 


(l  -  v2)Gc  a‘ 


^mn  sin  m^oc/a  sin  nny/b 

_ m2  +  n2(a2/b2) 

Qmn  sin  mioc/a  sin  rmy/b 
[m2  +  n2(a2/b2)J2 


For  the  case  of  a  uniform  load  intensity  q  =  Constant  and  for  the 
center  of  the  plate  (x  =  a/2,  y  =  h/2)  equation  (86)  hecomes 


Xb1+JL  ^f  (h  +  ■ 
VN  2  (l  -  v2)g.  a2 


sin  pht/2  sin  ng/2 
m[m2  +  n2(a2/h2)] 
sin  mjt/2  sin  ng/2 

infra2  +  n2(a2/b2)J2 


The  ratio  of  the  series  is  I.98,  when  a/b  =  1,  and  the  ratio  of  the 
series  is  1.11  when  a/b  =  l/2. 

For  the  case  of  a  concentrated  load  at  the  center  of  the  plate 
the  deflection  ratio  at  the  point  of  load  application  assumes  the  form 


1  *2 

1  +“ 


(h  +  t)t 


(1  -  v2)g, 


(sin  mrt/2  sin  mr/2)2 
m2  +  n2(a2/b2) 
(sin  mjr/2  sin  mt/2)2 

[m2  +  n2  (a2/h2  )J  2 


Now  it  is  easily  shown  that  the  numerator  series  in  equation  (88) 
does  not  converge  and  consequently  w/wjj  =  00  in  this  case.  A  more 
detailed  consideration  shows  that  in  any  plate  theory  which  takes 
transverse  shear  deformation  into  account  the  deflection  under  the  point 
of  application  of  a  concentrated  load  must  become  infinite  in  contrast 
to  what  happens  when  transverse  shear  deformation  is  not  taken  into 
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account.  This  difference,  of  course,  vanishes  as  soon  as  the  load 
intensity  he comes  finite,  and  then  the  theory  with  transverse  shear 
deformation  taken  into  account  is  more  accurate  than  the  theory  which 
does  not  take  into  account  this  effect . 

For  the  sake  of  numerical  illustration  take  again  the  square 
•plate  (a A  =  l)  with  uniform  load  distribution.  According  to  equation  (57) 
the  deflection  at  the  center  is  increased  because  of  transverse  shear  by 
the  factor 


JL  .  1  +  9.7  (h  +gtlk  (87a) 

wu  Gc  a^ 


Take  h  =  1.0  inch,  t  =  0.1  inch,  a  -  10  inches.  Ef/Gc  =200 
and  V  =  l/3.  Then,  according  to  equation  (87a),  v/wfl  =  1  +  2.3,  so 
that  in  this 'case  the  deflection  with  transverse  shear  is  more  than  three 
times  the  deflection  when  shear  deformation  in  the  core  is  neglected. 

Returning  now  to  equation  (84)  for  w  and  equation  (83)  for  <a 
and  substituting  these  two  equations  in  equations  (76)  and  (77)  in  order 
to  determine  the  changes  of  slope  and  after  slight  transfor¬ 

mations  there  results 


- cos  Xmi  sin  |iny 

(Xm2  + 

> 


Py 


+  M-n^) 


sin  cos  p.ny 


(8  9) 


Equations  (89)  are  remarkable  for  the  reason  that  they  are  not  affected 
by  transverse  shear  deformability.  According  to  equations  (73),  the 
same  is  then  true  of  the  bending  and  twisting  couples  %,  %,  and  Mjcy* 

It  is  not  easy  to  see  why,  in  this  statically  indeterminate  problem,  the 
magnitude  of  the  internal  forces,  as  well  as  that  of  the  deflections,  does 
not  depend  on  the  elastic  properties  of  the  core.  The  analysis,  however, 
shows  that  the  distributions  of  M*,  My,  and  M^,  and  therewith 
of  Qr/r  and  Qy,  remain  the  same  as  those  obtained  under  the  assumption 
that  Gc  =  00.  In  this  connection  the  following  remark  may  be  made. 
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Evidently  the  following  three  boundary  conditions,  w  =  %  =  0y  =  0 

along  the  edges  x  =  0,a,  have  been  satisfied.  In  order  that  the  last 
of  these  three  conditions  be  satisfied  there  are  necessarily  nonvanishing 
edge  values  of  the  twisting  couples  M^y.  The  same  is  true  in  the  theory 

without  transverse  shear  deformation,  where,  however,  no  alternative 
possibility  exists,  as  in  that  theory  only  the  boundary  condi¬ 
tions  w  =  =  0  are  relevant.  For  the  present  system  of  equations 

three  boundary  conditions  must  be  formulated  for  every  plate  edge.  Thus, 
it  is  possible  although  mathematically  complicated  to  solve  the  problem 
of  the  rectangular  simply  supported  plate  with  the  edge  condition  =  0 
replaced  by  the  condition  =0.  In  that  case,  which  will  not  be 

pursued  here,  there  evidently  will  be  a  distribution  of  internal  stresses 
which  is  modified  by  the  effect  of  transverse  shear  deformation 

Cylindrical  bending  of  plates.—  As  a  further  relatively  simple 
example  of  application  of  equations  (70),  (72),  and  (73)  problems  are 
considered  for  which 


d(  )/Sy  =  0 

d(  )/bx  =  d(  )/dx  =  (  )' 
Mjy  =  Qjr  =  m.y  =  =  0 


^  =  VM* 


(90) 


and  where  consequently  the  problem  reduces  to  the  following  system  of 
equations : 9 


Qj*  +  q  =  0 


—  Qj  +  —  0 


Or  =  (h  +  t)Gc(pI  +  w») 

(l  -  v2)Mr  =  D*Pxf 

^Note  that  in  order  to  obtain  the  problem  of  the  sandwich  beam  from 
equations  (90)  and  (91)  the  only  changes  which  are  necessary  amount  to 
setting  v  =  0  in  equations  (91) . 
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To  set  into  evidence  the  effect  of  finite  values  of  Gc  in 
equation  (91) ,  the  following  system  of  equations  is  deduced  from 
equation  (91): 


_  (IV)  ,  Dq" 

Dwv  '  =  q  +  mx‘ - - — 


(h  +  t)G( 


M*  =  -  Dwm  - 


Dq 


(h  +  t)G< 


%■  =  -  Dw’"  -  - — —  +  mx 
(h  +  t)Gc 


pT  =  -  w*  + 


(h  +  t)G( 


(92) 


(93) 

(94) 

(95) 


Solutions  to  the  following  problems  are  listed: 

(l)  Simply  supported  plate  of  span  Z  carrying  a 
load  q  =  qn  cos  rtx/Z.  Boundary  conditions:  w(±Z/2)  =  ^(±1/2)  -  0. 


w 


It2  Ef 

(h  +  t)t 

cos  nx/Z 

2  (l-V2)^ 

Z2 

U/Dk 

(96)10 


As  the  problem  is  statically  determinate  as  far  as  moment  and  force  are 
concerned  there  is  no  modification  of  M*-  and  Qx  due  to  the  finite 

value  of  Grc . 


^The  factor  in  brackets  may  again  be  written  in  the  form  1  +5.^0, 

with  3  =  (Ef/Gc)/[(h  +  t)t/Z2], 
reference  5. 


using  the  notation  suggested  in 
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(2)  Simply  supported  plate  of  snan  Z  carrying  a  uniform. 
load  q  =  qn. 


v  =  _  i]  _  i[i  +  SD  lf/xf  _ 

D16  ]2kp/2)  *]  (h  +  t)Gcl^\l^) 


"From  this  for  the  center  deflection. 


w(0)  =  Jfr  iOEfi  +.  24 _ l 

384  D  5  (l  - 


v2)g, 


It  is  seen  that  the  correction  factor  for  the  center  deflection  is  almost 
the  same  as  that  for  the  cosine  load  curve  (equation  (96)),  the  only 
difference  being  a  change  of  the  factor  n2/ 2  =  4.93  into  24/5  =  4.80, 
that  is,  a  reduction  of  the  shear  correction  factor  by  at  most  3  percent 
is  present.il 

(3)  Built-in  plate  of  span  Z  carrying  a  uniform  load  q  =  qg. 

The  boundary  conditions  are:  w(±Z/2)  =  PY(±z/2)  =  0  (and  ™ 
not  w’(±Z/2)  =0). 


_  <10* 


(h  +  t)GcZ2 


3  ($)  - 


From  this  there  follows  for  the  center  deflection. 


v(0)  = 


+  2i,  Sk+M 


384D I  (1  _  v2 )n  Z2 


(100) 


Comparison  of  equations  (100)  and  (98)  shows  that  for  the  built-in 
plate  the  effect  of  transverse  shear  deformation  is  very  much  more 


i:LNote  that  according  to  equation  (87a)  the  shear  correction  factor 
for  the  square  plate  of  width  a  =  Z  is  more  than  twice  as  large  as  the 
shear  correction  factor  for  the  plate  strip  of  width  Z. 
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pronounced  than  it  is  for  the  simply  supported  plate,  a  factor  24/5 
in  the  latter  case  being  replaced  by  a  factor  24  in  the  former  case. 12 
As  a  further  result  in  this  problem  of  the  built-in  plate,  by  putting 
equation  (99)  into  equation  (93) >  it  is  found  that  the  moment 
function  does  not  contain  any  terms  depending  on  the  effect  of 

transverse  shear  deformation.  This  again  is  somewhat  surprising  as 
in  this  case  it  is  not  possible  to  determine  the  moment  function  by 
statics  alone. 13 

Circular  plates;  rotational  symmetry.—  As  no  examples  of  solutions 
of  circular  sandwich— plate  problems  have  as  yet  been  published  and  as 
it  is  of  some  interest  to  determine  in  which  way  the  shear  correction 
factors  change  in  going  from  a  problem  for  the  plate  strip  to  the 
corresponding  problem  for  the  circular  plate,  the  equations  for 
ax i symmetrical  transverse  bending  of  circular  plates  are  briefly 
discussed. 

Polar  coordinates  r ,0  are  introduced  and  notation  which  is 
customary  in  plate  theory  is  used.  As  a  consequence  of  equa¬ 
tions  (70),  (72),  and  (73),  the  following  system  of  equations  is 
obtained: 


drQr 

—  +  r<* 


0 


dr  Mr 
dr 


Mg  -  rQp  +  rn^  =  0  J 


don 


^=(11  +  t)Gc(pr  +  dv/dr) 


(102) 


A  somewhat  similar  percentage  increase  must  take  place  in  going 
from  equation  (86)  for  the  rectangular  plate  with  all  four  edges  simply 
supported  to  a  formula  (which  has  not  yet  been  derived)  for  the  rectan¬ 
gular  plate  with  all  four  edges  built  in. 

!3as  a  problem  where  the  moment  distribution  is  in  fact  dependent 
on  the  effect  of  transverse  shear  there  may  be  mentioned  the  problem  of 
the  cylindrically  bent  plate  with  both  ends  built  in,  which  carries  a 
load  qj  =  q_jX  instead  of  the  load  q0  =  qQ.  This  problem  also  may  be 

solved  by  means  of  equations  (92)  to  (95)* 
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Mp  -  V%  =  D*  dPr/dr 
%  -  vMp  =  D*pr/r 


(103) 


According  to  equation  (79),  the  equation  for  the  deflection  w 
will  he 


I V^V^w  =  q 


1  1  drmp  _  y2q 

2  r  (h  +  t)Gc 


(104) 


where  =  (l/r)d[r  d(  )/dr]/dr.  Having  found  w  by  means  of 
equation  (104),  0r  may  he  determined  from 


0r 


dw  +  Or 

**  (h  +  t)Grc 


dw 

dr 


+ 


rq  dr 


(105) 


and  therewith  %  and  are  obtained  from  equations  (103). 

In  the  present  problem  it  seems  to  he  somewhat  more  convenient  to 
proceed  as  follows:  Combine  equations  (101)  and  (103)  to  obtain  as 
equation  for  the  change  of  slope  pr. 


D 


1  4_ 

r  dr 


<L  A 


r  dr\^  drA')_ 


=  1 


_  1 

r  dr 


(106) 


Having  Pr,  Mp  and  %  are  found  from  equation  (103)  and  Qr,  from  the 
second  of  equations  (101)^ 

«e  -  nft  Ur  +  V  j£t) 


(107) 
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(108) 

Finally ^  vith  this  value  of  w  is  found  “by  integrating  equation  (102) 


-I 


Pr  dr  + 


(h  +  t)G( 


i  drPA  +  Jmr  ^ 

?  dr  J  (h  +  t)G( 


(109) 


Deflection  of  circular  plate  vith  -built-in  edge.-  The  tending  is 

now  considered  of  a  plate  with  transverse  load  q  =  qn(r/a)n  and 
with  1%.  =  0.  First,  from  equation  (106), 


r  t  .  In8-3 


DPr  -  »X  S  +  02©  +  03  J  iog.  f  +  (n  +  —  -+-^ 


(C3 


Attention  is  restricted  to  complete  plates  with  no  concentrated 
load  at  the  center,  and  consequently  it  is  necessary  to  set  Cg  =  c^ 

in  equation  (110) .  This  gives 


Dpr  =  ci  |  + 


(n  +  4)(n  +  2) 


(110a) 


Putting  equation  (110a)  into  equation  (109),  there. results  for  the 
transverse  deflection  w 


D  -  = 


cl, 

rr\2  (  qn&3 

(Z) 

n+4 

4*  C)i 

2 

(n  +  4) (n  +  2)2 

w 

+  D  2c,  +  tna,3_/&^ 

a2(h  +  t)GcL  (n  +  2)2^  - 


(111) 
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Taking  the  case  of  a  plate  of  radius  a  with  "built-in  edge,  that 
is,  with  the  "boundary  conditions 


Pj.(a)  =  w(a)  =  0  (112) 

there  results 


Dpr 


<lna3 

(n  +  4) (n  +  2)2 


(113) 


and 


Dw  = 


_  3na 


(n  +  2)2 


n+4 


i2(h  +  t)Gf 


(n  +  4)2 


fer  -  3 


1  (r/a)2  -  1 

2  n  +  4 


(114) 


From  equation  (ll4)  there  follows  for  the  deflection  at  the  center  of 
the  plate 


Dw(0) 


Sna 


2(n  +  2)(n  +  4)^ 


1  + 


(n  +  4)2  'E-p  (h  +  t)t 

n  +  2  (l-v2)oc  a2 


(115) 


Consider  the  following  special  cases: 

(l)  Uniform  load  distribution  qn  =  From  equation  (115)^  it 

follows  for  the  ratio  of  deflection  with  and  without  transverse  shear 
deformation  that 


w(0)  =  n  o  Ef  (h  +  t)t 

[v(0j  _  m  “  (1  -  V2)GC  a2 

Lrc=co 


(116) 
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Equation  (ll6)  may  "be  compared  with  equation  (100)  for  the  deflection 
of  the  infinite  plate  strip  of  width  l  with  built— in  edges. 

Setting  I  «  2a,  it  is  seen  that,  while  the  transverse  shear  correction 
factor  for  the  strip  has  a  value  6,  the  corresponding  factor  for  the 
circular  plate  is  8.  This  is  consistent  with  the  earlier  comparison 
between  the  simply  supported  Btrip  and  the  simply  supported  square 
plate,  except  that  there  the  change  is  from  4.8  to  9*7. 

(2)  Linearly  increasing  load  distribution  q  =  qqr/a.  From 
equation  (115),  it  follows  that 


w(0) 

[w(0)]G 


(h  +  t)t 


3  (1  -  v2)o„  * 


(117) 


showing  that  the  correction  effect  iB  only  slightly  greater  than  in  the 
case  of  the  uniform  load  distribution. 

(3)  Load  increasing  linearly  from  edge  to  center,  q  =  qp  +  qq(r/a) 
/qi  =  -q0T!  From  equation  (115),  it  follows  by  superposition  that 


Dw(0)  = 


IQ**  ,  ,  o  Ef  (h  +  t)t 
2  X  2  X  Id  (l-v2)G,  a2 


2  x  3  x  25 


l  +^  — 

L  3  ci 


Ef  (h  +  t)t 
-v2)gc  ^ 


Wo)  =  43q°a  fi  +  22QO  Ef _ (h  +  t)t 

{>  32  X  150  387  (i_v2)gc  a2 


(118) 


Comparing  the  factor  3OOO/387  =  7.76  which  occurs  in  equation  (ll8) 
with  the  corresponding  factors  8  and  8.33  in  equations  (ll6)  and  (117) 
it  is  seen  that,  in  the  foregoing  three  problems  at  least,  there  is 
little  difference  between  the  transverse  shear  stress  correction  factors 
in  the  case  of  three  different  loading  conditions  for  the  circular, 
clamped-edge  plate.  The  fact  that  this  agreement  should  not  be  expected 
to  hold  generally  follows  again  by  considering  the  case  of  a  point 
load  at  the  center  of  the  plate,  for  which  the  shear  correction  factor 
would  again  be  infinite. 
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The  examples  of  this  section  should  "be  augmented  "by  the  solution 
for  the  circular  plate  of  radius  a,  which  carries  a  load  distributed 
uniformly  over  a  smaller  circle  which  is  concentric  with  the  boundary 
of  the  plate. 


Circular  Rings 

As  the  simplest  example  of  a  curved  sandwich  structure  there  are 
considered  in  this  section  stresses  and  deformations  of  circular  rings 
in  their  own  plane.  As  was  found  in  the  general  developments  of 
part  I  of  this  report,  in  a  curved  sandwich  structure  there  will  be 
the  effect  of  both  transverse  shear  and  normal  stress  deformation. 1^- 

There  are  set  for  the  relevant  coordinates  and  variables 


?1  “ 

ai  =  1 

El  -  a 

ui  =  V 

Pi  =  P 

£5 

11 

rH 

11 

Mj  -}  ss  M 

pi  =  p 

in^  =  m 

b(  ) /df  =  d(  ) /a  d 6  = 

(  )*/a 

H  -  X  =  |[(h  +  t)t/a2]  (Ef/Ec)  J 

The  equilibrium  equations  (25)  to  (28),  (31),  and  (3*0  reduce  to 
the  following  equations: 


*\ 

N’  +  Q  +  ap  =  0 
Q’  —  N  +  aq  =  0 


(120) 


l^The  effect  of  transverse  shear  stress  deformation  on  homogeneous 
circular  rings  has  "been  considered  hy  L.  Beskin  in  reference  6. 
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M'  —  aQ  +  am  =  0 
CT£m  =  8  ~  [M/(h  +  t)a] 


(121) 


The  stress— strain  relations  (56),  (59),  (60),  (6l),  and  (64)  to  (67) 
reduce  to  the  following  equations: 


(l  +  i  x)n  =  i  C[V  +  »  + 


(122) 


Q  =  (h  +  t)Gc|p  +  i(w*  -  v)J 


(123) 


(1  +  X)M  =  i  D*(P*  +  b/Ec) 


(124) 


The  load  terms  p,  q,  m,  and  s  are  given,  according  to 
equations  (19)  to  (22),  by 


(1  +  ^>+  (x-\* 


+  t 


PZ 


(1  +  -  (1  -  ^i)P; 


?  (125) 


h  +  t 


1 

8  =2 


A  h  +  t\  f 

i  h  +  t\ 

i1  +  2a  N  (, 

1  2a  fa 

^3 
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Ring  sector  acted  upon  by  end  bending  moments As  a  first  problem 
on  circular  rings,  which  illustrates  the  effect  of  transverse  normal 
stress  deformation,  there  is  taken  this  basic  case  for  which,  as  is 
known,  there  must  be  the  same  stress  distribution  at  all  sections 
9  =  Constant  of  the  ring. 

According  to  equations  (120)  and  (121), 


M  =  Mq 


(126) 


Equations  (122)  to  (l2h)  become 


v»  +  w  =  0  1 


P  +  (w*  —  v) /a  =  0 


(127) 


(1  +  =  D*P' /a 


(128) 


The  significant  result  of  this  consideration  is  contained  in  equation  (128), 
which  may  be  written  in  the  alternate  form 


M=  M, 

^  a  1  +  X. 


D* 


V 


1 1 


+  w 


W 


(h  +  t)t  Ef 
a2  E,, 


(129) 


Thus,  in  this  case  of  pure  bending  the  transverse  flexibility  of 
the  core  is  responsible  for  a  reduction  of  the  bending  stiffness 

factor  .D*  =  |  t(h  +  t)2Ef  which  is  obtained  exactly  when  Ec  =  0 

and  practically  when  Ec  is  of  the  same  order  of  magnitude  as  Ef . 
Equation  (129)  shows  that  the  reduction  of  D*  is  significant 
whenever  Ec  is  so  small  that  the  ratio  Ec/Ef  is  of  the  same  order 

of  magnitude  as  the  ratio  (h  +  t)t/a2. 
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As  a  numerical  example  take  the  following  values:  h  =  0.9  inch, 
t  =  0.05  inch,  a  =  20  inches,  and  Ef/Ec  =  1000,  for  which 


1  (h  +  t)t  Ef 

2  a2  Ec 


1  v  9j22  x  °‘°^ 

2  A  i+oo 


x  1000  =  0.0595 


indicating  a  reduction  in  tending  stiffness  of  about  6  percent. 

Changing  a  from  20  inches  to  10  inches  changes  the  effect  from 
6  percent  to  2k  percent.  Changing  Ef/E c  from  1000  to  2000  increases 
the  effect  from  6  percent  to  12  percent.  Altogether  it  may  he  said 
that  this  effect  is  of  noticeable  magnitude  for  some  geometrically 
reasonable  structures  when  the  modulus  ratio  Ef/Ec  is  of  the 
order  1000  or  more.  Assuming  aluminum  face  layers  with  Ef  =  10?  psi, 
this  means  that  Ec  ft?  lO1*’  psi,  which  is  well  within  the  range  of  some 
present>-day  core— layer  materials. 

Comparing  equation  (129)  with  the  earlier  formulas  for  the  effect 
of  transverse  shear  stress  deformation,  for  instance  with  equation  (ll6) 

in  which  a  represents  the  plate  radius  and  observing  that  Gc  ~  ^  Ec, 

it  is  seen  that  the  correction  terms  are  of  the  same  form,  the  difference 
being  an  appreciably  larger  numerical  factor  in  the  expression  representing 
the  shear  effect. 

Closed  circular  ring  acted  upon  by  uniform  radial  load.—  Having 
rotational  symmetry,  d/d@  =  0  and  v  =  0  =  0.  Also  set  p  =  m  =  0. 

The  remaining  equations  permit  the  determination  of  the  stresses  in 
the  face  and  core  layers  in  a  way  which  depends  on  the  extent  to  which 
the  load  is  applied  to  the  outer  (upper)  and  inner  (lower)  face  membranes. 
Equation  (12)  becomes 


(130) 


From  equations  (121),  it  follows  that 


=  s  “  [m/ ( h  +  t)a] 


(131) 
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The  stress— strain  relations  (122 )  to  (124)  give 

3£  =  _L(i  +  h\aQ  -  Ch  +  ^  -  Bill  +  2fc  _  (h  +  t)C*l  _  ag  ,  1  s 

a  C*^X  +  3 n  ^ca  "  C*[  +  3  12EC&2  J  "  C*  (132) 

and 


M  = 


D*  s 
1+1  aEc 


(133) 


A  closed  circular  ring  subjected  to  a  uniform  radial  load 
distribution  q  is  stressed  not  only  by  a  uniform  axial  force  N  =  ag, 
as  would  be  expected,  but  in  addition  is  stressed  by  a  uniform  bending 
moment  M,  the  magnitude  of  which  is  given  by  eguation  (133) .  The 
explanation  of  this  result  is  that  for  a  ring  with  relatively  soft 
core  the  circumferential  stress  distribution  depends  on  the  extent  to 
which  the  external  radial  load  is  applied  to  the  inner  and  outer  forces, 
respectively.  Roughly  speaking,  for  a  sufficiently  flexible  core  layer 
the  load  gu  goes  predominantly  into  the  outer  face  layer,  while  the 
load  qj  goes  predominantly  into  the  inner  face  layer. 

According  to  eguations  (46),  in  the  present  case  for  the  stresses 
in  the  two  face  layers. 


■  i  H  +  rh  M 


N,  =  i  N - 1 


h  +  t 


M 


(134) 


According  to  eguations  (130)  and  (133)  and  in  view  of  the  definitions 
of  D*  and  X,  this  may  be  written 


N, 


r  _  §./  ,  2  Is  \ 

u  ’al? + 1  +  y 


", '  f  (9 


,  2Xb  ) 
1  +  Xy 


(135) 
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Combining  next  equations  (131)  and  (133),  for  the  transverse 
normal  stress  in  the  core  layer,  the  following  expression  is  obtained: 


a£m  = 


s 

1  +  X 


(136) 


For  a  specific  example  assume  that  the  radial  load  is  applied 
entirely  to  the  inner  face  of  the  ring  so  that  qu  =  0  and,  according  to 

equation  (125), 


1  -  i1  - 
■  -  -  Ki  - 


(137) 


With  q  and  s  given  hy  equations  (137),  equations  (135)  and  (136) 
■become 


1  _  h  +  t\  aaz 

2  If  2a  Jl  +  X 

if,  h  +  tV1  +  2X>a9z 
2\  2a  )  1  +  X 


(138) 


°5m  = 


h  +  t\  gl 
2a  Jl  +  X 


(139) 


It  is  seen  that  the  flexibility  of  the  core  layer  increases  the  circumfer¬ 
ential  stress  in  the  loaded  face  layer  in  the  ratio  (l  +  2\)/(l  +  X)  and 
decreases  it  in  the  unloaded  face  layer  in  the  ratio  l/(l  +  X) ,  where  X 
is  defined  by  equation  (119),  compared  with  the  equal  values  of  these 
stresses  when  Ec  =  00. 

Considering  once  more  the  numerical  data  under  the  section  entitled 
"King  sector  acted  upon  by  end  bending  moments,"  it  is  found,  for 
instance,  that  the  stress  in  the  inner  face  layer  may  be  about  6  or  12, 
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or  24  percent  higher  than  the  corresponding  stress  calculated  without 
taking  into  account  the  transverse  flexibility  of  the  core  layer. 

Ring  sector  acted  upon  hy  radial  loads  qu  and  qjj  uniform  in 
circumferential  direction  and  with  vanishing  resultant  q.-  Again  it  is 
assumed  that  d(  ) /d©  =  0,  m  =  p  =  0  and  now  in  addition  that  q  =  0, 
so  that,  according  to  equation  (125),  the  only  nonvanishing  load 
term  is  s.  Further  it  is  assumed  that  the  ends  8  =  ±  a  of  the  ring 
sector  are  free  of  stress,  that  is,  N(±a)  =  Q(la)  =  M(ia)  =  0.  The 
ordinary  theory  of  circular  rings  would  then  indicate  the  absence  of 
deformations  in  the  entire  ring.  In  the  present  case  there  is  found 
a  type  of  deformation  peculiar  to  the  sandwich  ring,  which  may  perhaps 
be  compared  to  the  action  of  a  Bourdon  gage. 

Solving  first  equations  (120)  and  (121)  and  satisfying  the  end 
conditions  of  the  ring  sector. 


N  =  Q  =  M  =  0 

=  8 


(l4o) 


The  stress— strain  relations  (122)  to  (124)  are  then 


vT  +  w  =  0 


aP+w*  —  v  =  0  > 


(141) 


p»  =  — s/Ec 


Assuming  s  independent  of  0,  from  equation  (l4l)  there  is  obtained  by 
integration,  with  constants  of  integration  Aj,  A^,  and  A^, 


p  =  -  A  0  +  A! 

v  =  — a  8  +  Aja  +  Ar>  cos  8  +  A3  sin  8 

c 

w  =  a  -S-  +  Ap  sin  8  —  Ao  cos  8 


> 


(142) 
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As  a  specific  example  consider  a  complete  ring,  slitted  radially 
at  the  section  0  =  it,  so  that  a  =  rr.  Prescribe  furthermore  the 
symmetry  conditions  P(0)  =  v(0)  =  w(0)  =  0.  Under  these  conditions 
there  is  obtained  from  equation  (l42) 


EqP  =  — 36 


Ecv  =  -as(0  —  sin  0) 


> 


(143) 


Ecw  =  as(l  —  cos  0) 


From  equations  (143),  it  follows  that  the  radial  slit,  which  is 
of  zero  width  before  the  loads  qu  and  qj  are  applied,  opens  under 
the  action  of  the  loads  to  a  width  given  by 


v(-jt)  -  ▼(*)  =  2*a  ^  -  2*a(l  +  (l^) 


For  a  numerical  example  take  a  =  10  inches,  h  =  1  inch,  t  =  0.05  inch, 
Fc  =  10,000  psi,  and  qu  =  20  psi,  and  obtain 

v(— it)  —  v(jt)  =  0.132  inch  (1^5) 


hare  “been  discussed 
simply  the  effect  of 
of  curved  sandwich 
effect  of  transverse 

Bending  of  semicircular  ring  by  end  shear  forces.—  Now  a  problem  is 
considered  in  which  both  the  values  of  Ec  and  Gc  affect  the  result 
of  the  analysis.  In  the  equilibrium  equations  (120)  and  (121)  all 
external  load  terms  are  set  equal  to  zero  and  then,  by  integration  and 
from  the  boundary  conditions,  that  is,  from 


The  foregoing  three  examples  of  ring  analysis 
in  some  detail,  because  they  illustrate  relatively 
transverse  normal  stress  deformation  in  the  theory 
structures,  without  involving  at  the  same  time  the 
shear  stress  deformation. 


N(±  f)  =  m(±  f )  =  0 

Q(k|)  =  iQo 


(146) 
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the  following  expressions  for  N,  M,  and 


Q  are  obtained: 


Q  =  Qo  sin  9 
N  =  Qq  cob  9 
M  =  -aQo  cos  9 


The  stress— strain  relations  (122)  to  (124)  become 


(l  +  X/3)Qo  cos 
Qq  sin  9  =  (h  + 
— (l  +  X)Qoa  cos 


9  =  (C*/a) (v*  +  w) 
t)Gc|p  +  (w*  -  v)/a] 
9  =  (D*/a)P* 


Integration  of  the  last  of  equations  (l48)  gives 

D*p  =  -a2(l  +  X)Q Q  sin  9 


(147) 


(148) 


(149) 


where  a  constant  of  integration  has  been  eliminated  by  means  of  the 
symmetry  condition  0(0)  =  0.  Substituting  equation  (149)  in  the  second 
of  equations  (l48) , 


!(»•  -  v)  -  Qo  Bin  x)] 

a2  fn  .  1  (h  +  t)t/Ef  Ef\l 

S'L1  +  5  (ft  +  ftjl 


=  Qq  sin  9 


(150) 


Simultaneous  solution  of  equation  (150)  and  the  first  of 
equations  (l48)  for  v  and  w  gives  as  general  expressions  for  v 
and 


v  =  A 0  cos  9  +  Aj  sin  6  +  A^  cos  9 


t 


W  =  A 9  sin  9  — 


+  B)cos  9  +  Ag  sin  9 


(151) 
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where  Ap  and  A2  are  arbitrary  constants  of  integration  and  A 
and  B  are  found  to  he 


,  1  (h  +  t)t/Sf 

2D*  2  a2  \EC  G  CJ 


1  (h  +  t)t/E-f  Ef\ 

2  &  +  G~c) 


(152) 


As  further  conditions,  it  is  prescribed  that  v(0)  =  =  0,  which 

makes  Ag  =  Aj  =  0  in  equation  (151).  There  remains 


v  —  AO  cos  6 


r 


w  =  AG  sin  0  —  B  cos  6 


(153) 


Of  particular  interest  are  the  values  of  w(jt/2)  and  w(0),  the  first 
of  these  giving  the  radial  deflection  of  the  point  of  load  application^ 
the  second  giving  the  change  of  radius  at  right  angles  to  the  applied 
load.  It  is  found  that 


(15^) 


(155) 


Equations  (154)  and  (155)  contain  the  interesting  result  that,  for 
this  problem,  transverse  shear  and  transverse  normal  stress  affect  the 
outcome  formally  in  nearly  the  same  way.  If  the  generally  unimportant 
terms  with  D*/a2C*  are  omitted,  which  amounts  to  the  usual  assumption 
of  circumferential  inextensibility  of  the  ring,  then  the  effects  of 
finite  Ec  and  Gc  occur  in  exactly  the  same  way. 
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For  a  numerical  example  take  h  =0.9  inch,  t  =  0.05  inch, 
a  =  20  inches,  Ef/Ec  =  1000,  and  Ep/Gc  =  2000.  This  gives 


1  (h  +  t)t  =  1 

2  a2  1^,800 

D*  =  1  fh  +  t)2  _  1 

a2C*  I  a2  1770 

,  =  1  (h  +  t)t  %  _  1 

2  a2  Ec  lO 

1  fh  +  t)t  %  =  2 

2  a2  Gc  16.8 

The  factors  in  "brackets  in  equations  (15^-)  and  (155)  "become 
1  +  16.8  +  16.8  +  1770 (X  +  3  X  1 6.8)  ~  1,18 

and 

1  +  +  -2— - (i  + - =  1.18 

16.8  16.8  1770  V  3  x  16.8/ 


Thus,  in  the  present  example  the  flexibility  of  the  core  is  responsible 
for  an  l8-percent  increase  of  deflection-load  ratio,  and  of  this  12  percent 
is  due  to  transverse  shearing  and  6  percent  to  transverse  normal  stress. 
Compared  with  these  two  effects  the  effect  of  circumferential  extensibility 
of  the  composite  ring  is  seen  to  be  negligible.  As  a  further  numerical 
illustration,  it  is  noted  that  reducing  the  ring  radius  a  from  20  inches 
to  10  inches,  with  all  other  data  unchanged,  changes  the  l8-percent 
correction  to  a  72—percent  correction. 

Bending  of  complete  circular  ring  under  action  of  two  concentrated 
radial^ forces  at  6  =  ±n/2 .—  The  solution  of  this  problem  may  be~ "obtained 
by  superposition  of  the  solutions  for  the  semicircular  ring  under  the 
action  of  end  shear  forces  Qq  (equations  (l46)  to  (155))  and  under  the 
action  of  end  loading  moments  Mq  (equations  (126)  to  (129)). 
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The  first  step  consists  in  determining  Mq  in  terms  of  Qq  such 
that  the  sum  of  the  P’s  from  equations  (129)  and  (1^9)  assumes  the 
value  zero  for  9  =  it/2,  that  is,  the  value  of  the  superimposed  "bending 
moment  at  9  =  it/2  must  make  the  tangent  to  the  deflected  ring  at 
this  point  horizontal.  Combining  equations  (129)  and  (1^9)  in  this 
manner,  there  is  obtained 


JT  1  +  1 

2  D* 


aMo 


Is-2  -  0 


or 


Mq  =  (2/ir)aQo 


(156) 


It  may  "be  noted  that  equation  (156)  is  a  further  case  of  a  statically 
indeterminate  problem  where  transverse  shear  and  normal  stress  flexibility 
do  not  affect  the  internal  force  and  moment  distribution  but  affect  only 
the  state  of  deformation  of  the  structure. 

Further,  the  radial  deflections  w(jt/2)  and  v(0)  due  to  the 
action  of  Mo  are  calculated,  in  order  to  combine  them  with  equations  (152) 
and  (153).  Integrating  equations  (129)  and  (127)  with  the  boundary 
conditions  v(0)  =  v(it/2)  =  0,  there  is  obtained  for  the  displacements 
due  to  Mo, 


D*w  =  — (l  +  l)MQa2^l  —  ^  cos  9^ 
D*v  =  (1  +  l)MQa2^0  -  £  sin  e'j 


y 


(157) 


and^  in  particular^ 


D*w(0)  =  (l  +  l)Moa2^  -  l) 
D*w(|)  =  -(1  +  l)Moa2 


(158) 
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Combining  equations  (158)  with  equations  (154)  and  ( 155)  and  taking  Mq 
from  equation  (156),  there  follows  for  the  resultant  displacements 


where  Xn  = 


(h  +  t)t  Ef 
2a2  Gc 


has  been  put  as  a  further  abbreviation. 


Equations  (159)  may  "be  written  in  thei  alternate  form 


(160) 


(161) 


inextensible,  which  amounts  to  putting  D*/a^C*  =  0  in  equations  (l60) 
and  (l6l),  then  equations  (l6o)  and  (l6l)  ' reduce  to  well-known  results 
of  circular-ring  analysis. 

Comparing  equations  (l60)  and  (l6l)  for  the  closed  circular  ring 
with  equations  (154)  and  (155)  for  the  open  semicircular  ring,  it  is 
noteworthy  that  for  the  semicircular  ring  X  anri  Xq  occur  with  equal 
weight,  while  for  the  closed  circular  ring  the  influence  of  Xq  is 
considerably  greater  than  the  influence  of  X.  Thus,  for  the  closed 
circular  ring  the  effect  of  transverse  shear  deformation  is  much  more 
important  than  the  effect  of  transverse  normal  stress  deformation,  while 
for  the  open  semicircular  ring  both  effects  occur  in  a  much  more  nearly 
equally  important  way. 

For  a  numerical  example  of  the  use  of  equations  (l60)  and  (l6l) 
take  again  the  values  for  the  numerical  example  given  in  the  section 
entitled  "Bending  of  semicircular  ring  by  end  shear  forces."  This 
gives  for  the  expressions  in  braces 


w(*/2)  =  0.149  | 


1  +  x  +  5.29[xa  +  Jl-(1  +  1)] 


w(0)  =  -0.137  -j 


1  +  x  +  3.65  k 


When  X  =  Xq  =  0  and  when  the  composite  ring  is  assumed  axially 
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!  +  +  2  X  1*22  +  =  l  69 

X  +  16.8  16.8  +  1770  y 


and 


1  + 


16.8 


2  X  3.65  .  3.65 


16.8 


1770 


=  1.50 


Thus,  while  the  effect  of  transverse  stress  deformation  for  the  open 
circular  ring  amounted  to  18  percent,  the  corresponding  corrections 
for  the  closed  ring  are  69  and  50  percent,  respectively. 

The  next  step  in  the  analysis  of  sandwich— type  circular  rings 
would  "be  the  general  solution  of  the  system  of  equations  (120) 
to  (124)  for  arbitrary  load  distributions.  This,  evidently,  is 
possible  and  further  specific  examples  of  interest  might  be  analyzed 
on  the  basis  of  the  general  solution.  Such  extension  of  the  work  of 
this  section  is,  however,  left  for  future  considerations. 


Circular  Cylindrical  Shells 

In  this  section  the  general  system  of  equations  of  part  I  of  this 
report  is  restricted  to  the  equations  of  the  theory  of  circular  cylindrical 
shells.  The  treatment  of  sandwich— type  shells  of  this  kind  is  shown  to 
be  not  appreciably  more  difficult  than  the  analysis  without  the  effect 
of  transverse  shear  and  normal  stress. 

As  specific  examples  some  problems  of  rotationally  symmetric 
deformations  are  treated.  In  particular  the  influence  coefficients 
are  obtained  for  a  semi— infinite  shell  acted  upon  by  bending  moments 
and  transverse  forces  at  one  end  of  the  semi— inf inite  shell.  With  these 
influence  coefficients  an  explicit  solution  is  obtained  for  the  problem 
of  the  infinite  circular  cylindrical  shell  acted  upon  by  a  pressure 
band  of  zero  width. 


NACA  TN  No.  1832 


55 


In  the  general  equations  of  the  problem  there  are  set  for  the 
relevant  coordinates  and  variables. 


lx  =  aS 

i2  =  x 

=  a 

8 

11 

CVl 

Ph 

II 

1— 1 

ro3 

ro 

11 

Qi  =  Qs 

02  =  Qx 

M11  =  % 

£ 

II 

OJ 

U]_  =  V 

u2  =  u 

Pi  =  Ps 

P2  =  Px 

mi  =  “S 

II 

S' 

Pi  =  PS 

P2  =  Px 

ai 

=  02-1 

(162) 


W12  =  N21  =  WxS 


%2  =  J 

The  equilibrium  differential  equations  (25)  to  (28),  (3l),  and  (34) 
■become 


1  <^X6 

8x  a  ^9  x 


8n„ 


xS  +  1  + 


8x  a  8s 


Pe  -  0 


+  1  -  2s  +  4  =  0 

8x  a  80  a 


(163) 


1  ^s 

cht  a  <30 

dMgS  ,  1  ^8 
8x  a  8s 


—  %c  +  “be  -  0 

-  Qg  +  m0  =  0 


(164) 
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°U  =  8 


Mg 

(h  +  t)a 


(165) 


The  stress— strain  relations  (56),  (59) ,  ( 60) (6l),  and  (64)  to  (67) 
■become,  with  X2  =  ^12  =  0,  X^  =  ^|(h  +  t)t/a^J(lf/Ec)  s  X, 


(h  +  t)q 
12aEc 


Nz  -  VNe  = 

2(1  +  V)NZ0  = 


be) 


> 


(1 66) 


^  -  (h  +  t)Gc(p0  +  J  §|  -'?) 

Qj.  =  (h  +  t)Gc^pz  + 

D*(i  Hje  +  Jj) 

^  ^  “  Hs1) 

2(1  ♦  v)m*  -  B^|a  +  J  I*) 


(167) 


(168) 


When  Gc  =  Ec  =  »  (and  therewith  X  =  0)  equations  (163),  (l64), 
(l66),  (167),  and  (l68)  reduce  to  the  known  system  of  equations  in  which 
deformations  due  to  transverse  stresses  are  neglected.  The  solution  of 
the  present  system  of  equations  is  not  essentially  more  difficult  than 
the  solution  of  the  system  with  Gc  =  Ec  =  00.  In  particular  also  here 
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there  may  he  obtained  a  trigonometric  double-series  solution,  as  a 
generalization  of  Navier's  solution  for  the  flat  plate  (references  7 
and  8) . 

For  this  trigonometric  double-series  solution  there  is  set. 


q  =  ^  sin  m©  sin  nx/l 

Pq  =  ^  ^  Pemn  cos  m©  sin  nx/l 

Pi  =  ^  Pxmn  sin  cos  D*/l 

%:  =  ^  ^  111,^,  sin  m©  cos  nx/l 

^  m^jnn  cos  m ©  sin  nx/l 

s  =  Smn  s in  m©  sin  nx/l 


(169) 


w  =  ^  Vnm  sin  m©  sin  nx/l 

v  =  ^  Vrrm  cos  m©  sin  nx/l 

U  =  ^  Uj,^  sin  m©  cos  nx/l 
P*  =  ^  Pxmn  sin  m©  cos  nx/l 

P©  =  ^  ^  P©mn  008  tt©  sin  nx/l 


(170) 
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«r-  II  sin  m 6  cos  nx/Z 

Q@  =  ^7  %mn  cos  m0  sin  raA 

(Nx,  Ue)  =  y~  J~  (Wxmn,  N0mn)  sin  m<?  sin  nx/Z 

Nx0  =  II  Nx0mn  cos  m0  cos  nx/Z 

(*4c>  Me)  -  ^  Mgjnn)  sin  m0  sin 

MjcO  =  ^  ^  Mxeinn  cos  m0  cos 


(171) 


When  equations  (169)  to  (171)  are  substituted  in  equations  (163) 
to  (l68)  there  remains  for  every  value  of  m  and  n  a  system  of 
13  simultaneous  equations  for  the  13  Fourier  coefficients  which  occur 
in  equations*  (170)  and  (l7l). 

A  system  of  only  five  simultaneous  equations  for  the  five  Fourier 
coefficients  in  equation  (170)  is  obtained  if  first  equations  (163) 
and  (l64)  are  reduced  to  five  equations  for  the  five  unknowns  w_,  vy  u, 
and  by  means  of  equations  (l66)  to  (l68)  . 

For  the  present s  the  task  is  not  carried  out  of  obtaining  the 
deformation  and  internal  stress  Fourier  coefficients  of  equations  (170) 
and  (171)  in  terms  of  the  Fourier  coefficients  of  the  load  terms  in 
equation  (169).  Instead^  the  axisymmetrical  case_,  to  which 
equations  (169)  to  (171)  reduce  when  sin  m 0  and  cos  m 0  are 
interchanged  throughout ^  and  then  only  the  terms  for  m  =  0  are  taken^, 
is  treated  separately. 

Axisymmetrical  deformation  of  circular  cylindrical  shell.—  In 
equations  (163)  to  (l68)  set 


S(  )  fdG  =  0 
&(  )/Sx  =  (  )' 

WX0  =  %  =  %x  =  0 
v  =  P0  =  0 

m0  =  p0  =  0 


(172) 
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and  then  the  following  system  of  equations  has  to  he  dealt  with: 


+  Px  =  0 


V  -  (%/a)  +  q  =  0 


(173) 


Mjc*  ~  Qx  +  m*  =  0 
c^m  =  s  -  %/(h  +  t)a 


(174) 


(1  * 

Nz  -  VN0 


-  VNy  =  C* 


(h  +  t)ql 


12  aE, 


=  C*u' 


(175) 


Or  -  (h  +  t)Gc(Px  +  w») 


(176) 


(1  +  X)t%  -  VM^  =  D*s/aEc 
^  -  V%  = 


(177) 


The  system  of  equations  (173)  to  (177)  may  he  reduced  to  two 
simultaneous  equations  for  Px  and  0*.,  as  follows:  First, 

express  in  terms  of  Px  hy  means  of  equation  (175)  and  substitute 

the  result  in  equation  (174).  From  the  first  of  equations  (177),  it 
follows  that 


D*s 


(178) 


(l  +  X)aE0 


6o 
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and  this,  introduced  into  the  second  of  equations  (177),  gives 


M,c  = 


(1  +  X)D» 
1  +  X  -  V2 


V 


VD*s _ 

aEc(l  +  X  -  V2) 


(179) 


Equation  (179)  is  introduced  into  the  first  of  equations  (17*0  and, 
restricting  attention  to  shells  of  uniform  section  properties,  there 
is  obtained 


(1  +  x)d* 

1  +  X  —  V2 


-  Qjc  =  -%  - 


vp*s _ 

aEc(l  +  X  -  V2) 


(180)* 


To  obtain  the  second  of  these  equations,  first,  introduce  into 
equation  (176)  the  value  of  w’  which  follows  from  equation  (175), 
giving 


(h  +  t)Gc 


C*(h  +  t)q* 
12aEc 


(181) 


In  equation  (l8l),  N0*  and  Nx’  are  taken  from  equation  (173)  and, 
after  slight  transformations,  there  is  obtained 


C* 


(h  +  t)Gc 


+  0* 


(182)* 


Comparing  equations  (180)  and  (182)  with  the  corresponding 

equations  without  the  effect  of  transverse  shear  and  normal  stress 

deformation,  it  is  seen  that  the  effect  of  transverse  normal  stress, 

which  is  represented  by  X,  merely  somewhat  modifies  some  of  the 

coefficients  of  the  left  sides  of  the  corresponding  system  of  equations 

with  E„  =  00.  In  contrast  to  this,  the  effect  of  finite  Gc  is  to 
c 

introduce  a  new  term  into  the  left  Bides  of  these  equations.  This  new 
term  may  be  of  appreciable  importance,  as  will  be  shown. 

Having  solved  equations  (l8o)  and  (182),  and  Mg,  are  obtained 
from  equations  (179)  and  (178),  respectively;  N0  follows  from 
equation  (173)  in  the  form 


Ne  =  a(Qx*  +  q) 


(183) 
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and  w  follows  from  equation  (175)  in  the  form 


w  =  (a/C*) 


(l  +  X/^aQ^'  +  aq  +  v  ^  pz  dx 


The  following  examples  illustrate  the  use  of  equations 

Infinite  circular  cylindrical  shell  with  periodic  load 
In  specialization  of  equations  (169)  to  (171) ,  set 


q  =  q„  sin  \a 


8=8^  sin  nx 


Px  =  pztl  cos  mi 


%  =  COS  |ix 


w  =  w^  sin  pa 
u  =  u^  cos  mi  > 

Pz  =  pz^  cos  mi 

=  V  cos  ** 

Nz  =  NZ(i  sin  mi 

Ne  =  n@m  sin  [ 
Mjc  =  Mxm  sin  MI 
M©  =  *%M  sin  mi 


6l 

(184) 

(178)  to  (184). 
distribution 

(185) 


(186) 


(187) 
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By  introducing  equations  (185)  to  (187)  into  equations  (l8o) 
nnfl  (182) ,  two  simultaneous  equations  are  obtained  for  the  ampli¬ 
tudes  and  0X(i,  as  follows: 


lVx^v2  ^  +  =  +  1 


+  x  -  v2  ^c 


VD* 


-2 

3  _  a_ 

pxn  c* 


H) 


ii2  + 


2tE. 


a2(h  +  t)G, 


2 

V  =  -  ^(^)i  +  vPx4/a 


^  (188) 


To  simplify  the  further  discussion,  hy  setting  in  equation  (l88) 
=  pX|I  =  0,  there  is  obtained  for  ^xn  and 


_  1  +  1  -  V2  V 

(1  +  \)D*  ^3 


K 


The  quantity  K  is  given  by 


(189) 


K  = 


1  +  X  +  Q_  Z2t  % 
3  it2  (h  +  t)a2  Gc 


1  +  X  -  V2  kj} 

(1  +  X)jA  (h  +  t)2a2 


(190) 


where  use  has  heen  made  of  the  relation  p  =  jr/2.  In  equation  (190) 
the  term  X/3  will  usually  he  of  little  importance.  The  other  two 
variable  terms  represent  the  effect  of  transverse  shear  deformation 
and  of  shell  curvature,  respectively.  When  the  radius  a  is  so  large 

that  Z^/(h  +  t)^a^  «  1,  the  shell  behaves  under  the  action  of  the 
given  load  essentially  as  a  plate  strip.  The  effect  of  transverse 
shear  is  important  as  soon  as  the  term  (2/jt2)  ( Z2/a2)  (t/(h  +  t))(Ef/G-c) 
is  not  small  compared  with  1. 

Before  evaluating  a  numerical  example  the  following  further 
formulas  which  are  readily  obtained  from  equations  (179)*  (183)* 
and  (l84)  are  listed: 


MCA  TN  No.  1832 


63 


o 

a^CLn 

T J  —  -  ^ 

n  c* 


[1  -  (1  +  i)K 


(191) 


Equations  (191)  show  that  in  this  problem  not  only  is  the  deflection 
increased  because  of  the  effect  of  transverse  shear ,  and  with  that 
the  hoop  stress  resultant  Ng^  but  now  also  an  effect  is  found  on 

the  bending-moment  distribution  in  the  opposite  sense.  The  effect 

of  transverse  shear  is  to  reduce  the  magnitude  of  the  bending  moments 
in  the  shell.  This  result  is  in  contrast  to  what  was  found  for  the 
examples  which  were  worked  out  in  the  sections  on  plate  analysis  and 
circular  ring  analysis  and  is  therefore  of  particular  significance. 

Equation  (191)  for  w^  may  be  compared  with  the  corresponding 
expression  for  a  simply  supported  plate  strip  of  width  l,  with 
sinusoidal  load.  The  result  for  this  case  must  follow  from  equation  (191) 
in  the  limit  a-^-oo  and  agree  with  equation  (96) }  which  was  previously 
obtained.  To  compare  the  last  of  equations  (191)  with  equation  (96), 
the  last  of  equations  (191)  is  written  in  the  form 


1  + 


1  +  1  -  V2  n\k 
(1  +  X)D*  W 

1+  X  -  V2(l\k  C* 
(1  +  X)D*  W  a2 


1  +  uf  -1  +  x  1  (h  +  t)t  Ef 

ill  1  +  x  -  y2  2  z2  Q-c 

x  +  ( 1\ 2  2t  % 

n  3  W  h  +  t&, 

. 

1  +  ^  v  /2N  C* 

(1  +  \)D*  W  a2 


in  (192) 


Equation  (192 )  reduces  to  the  equivalent  of  equation  (96)  if  in 
it  a  — 3p>  oo# 

From  a  comparison  of  equations  (192)  and  (96),  it  is  further 
concluded  that  the  correction  due  to  transverse  shear  is  greatest 
in  this  case  when  a  =  oo5  so  that^  in  this  case^  the  curvature  of 
the  shell  tends  to  reduce  the  additional  shear  deformation,  below 
the  value  obtained  for  the  simply  supported  plate  strip.  ’ 
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For  a  numerical  example  first  take  h  =  1  inch,  t  =  0.05  inch, 
a  =  10  inches,  1-20  inches,  Ef/Gc  =  200,  Ef/Ec  =  100,  V  =  1/3, 

and  X  =  |  100  =  °-°25 .  The  factor  X  of  equation  (190) 

■becomes 


r  fi  n  nnn  2  400  X  0.05  „  „A  4(1-  0_.09l  160x000 
K  v  5  100  X  1.05  x  100  +  *4  100  X  1.: 


=  (1  +  0.008  +  3.86  +  54.5)  =  0.01685 


while  without  transverse  shear  and  normal  stress  deformation 


(K)Gc=Ec=c  =  (1  +  54. 5 )-1  =  0.0180 

The  correction  in  this  case  amounts  to  about  6  percent. 

Changing  the  moduli  ratio  to  Ef /G-q  =  2000,  Ef/Ec  =  1000, 


X  =  (1  +  0.08  +  38.6  +  54.5)  =0.0106 


instead  of  X  =  O.O1685.  The  correction  in  this  case  amounts  to 

0.0180  -  0.0106  x  100  M  ye  percent .  Thus  again  there  is  a  case  where 
0.0106 

omission  of  the  effect  of  transverse  shear  deformation  would  give 
results  which  could  not  he  used.  However,  it  is  noted  that  the  effect 
of  transverse  normal  stress  deformation  is  quite  small  and  may  here 
safely  he  neglected. 

If  the  foregoing  values  of  K  are  introduced  into  equations  (191) , 
it  is  seen  that  the  percentage  corrections  apply  to  the  bending-moment 
value  directly  hut  that  for  hoop  tension  and  radial  deflection  the 
corrections  are  very  small  indeed.  In  fact,  in  order  that  there  he 
appreciable  corrections  due  to  transverse  shear  on  hoop  tension  and 
radial  deflection,  it  is  necessary  that  the  half  wave  length  of  the 
sinusoidal  load  q  he  so  small  that  X  is  at  least  of  magnitude  0.25 


or  more. 
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A  case  of  approximately  this  kind  is  obtained  if  the  half  ware 
length  l  is  changed  from  20  inches  to  10  inches  and  the  moduli  ratios 
are  again  taken  as  Ef/Ec  =  100,  Ef/Gc  =  200.  Then, 


K  =  (1  +  0.008  +  O.965  +  3.4l)-1  =  O.1865 


whereas 


(k)Ec=Gc=00  =  (1  +  3  Mr1  =  0.227 

The  percentage  change  of  K  and  therewith  of  Mjr  is  slightly  more 
than  19.  The  percentage  change  of  %  and  w  is  about 

The  foregoing  numerical  examples  show  that  the  effect  of  transverse 
shear  may  be  significant  in  cylindrical  sandwich— shell  analysis  and  that 
moreover  its  magnitude  will  not  in  general  be  predictable  by  the  analysis 
of  an  equivalent  flat— plate  or  straight— beam  problem. 

For  the  infinite  circular  cylindrical  shell  with 
load  q  =  q^  cos  |jx  the  essential  results  are  given  by  equations  (190) 

and  (191).  These  results  may  be  extended  directly  to  the  loading 
condition 


1  = 


^  qn  cos 


M-nx 


> 


Mn  =  n*A  J 


(193) 


By  superposition,  from  equation  (l9l)  the  following  formulas  are 
obtained: 


_  ^ 

Mjc  =  Y_  COS  [lnx 

%  =  a  9n(1  -  Kn)  cos  hnx  > 

w  =  (a2/C*)  ^  qn jjL  -  (l  +  X/3)KnJ  cos  i^x 


(19^) 
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The  values  of  Kn  are  obtained,  from  the  formula 


K 


n 


,  ,  X  .  2  22t  %  4(1  +  X  -  V£)  2^ 

3  n2^2  (h  +  t)a2  Gc  (1  +  (h  +  t)2a2 


(195) 


Having  the  solution  for  the  infinite  shell  with  periodic  load 
distribution,  it  will  he  only  necessary  to  add  to  this  the  general 
solution  of  the  differential  equations  without  external  load  terms, 
in  order  to  obtain  the  complete  solution  for  any  edge  condition  of 
the  axisymmetrically  stressed  circular  cylindrical  shell  of  finite  length. 
This  additional  solution  will  now  he  obtained. 

Finite  circular  cylindrical  shell  acted  upon  by  edge  moments  and 
forces .—  To  solve  equations  (l8o)  and  (182)  with  right-hand  sides  equal 
to  zero,  equation  (182)  is  differentiated  twice  and  Px' '  is 
substituted  from  equation  (l80) .  This  gives 


a£/l  +  A 
C*V  3 


Qx 


IF 


1 

(h  +  t)Gc 


Q, 


I  t 


+ 


1  +  x  -  v2 

(1  +  x)d* 


or 


-  Smj2^’  *  +  4m21|Qx;  =  0 


(196) 


where 


mi 


V 


'i  c* 

(r  +  X/3)(h  +  t)Gc 


1  I  t  Ef 

a  \j (1  +  X/3)(h  +  t)  Gc 


(197a) 


JHq  = 


1  + 


x-v2 


a2D*  (1  +  X)(l  +  X/3) 


1  ,4 1  1  +  X  -  V2 

\J  (h  +  tjaV  +  +  V3) 


(197b) 
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Hie  auxiliary  equation  corresponding  to  equation  (19 6)  is 


r^  -  2m12r2  +  kva^  *  0  (198a) 

or 

r2  =  m-f  ±  (198b) 


The  solution  of  equation  (196)  occurs  in  two  different  forms,  depending 
on  whether  r2  of  equation  (198b)  is  real  or  not.  According  to 
equations  (197)  and  (198b),  r2  is  complex  as  long  as 


m^*  <  kmq 


or 


(199a) 


i_r _ tEf _ ] 

a^  _(h  +  t)(l  +  X/3)GC 


2 


< 


M]l  -  v2/d  +  \)] 

(h  +  t)2a2(l  +  X/3) 


To  clarify  this  condition,  neglect  X  (which  is  of  very  little  importance 
here)  and  equation  (199a)  then  becomes 


2a 

t 


When  equation  (199)  holds,  a  quantity  k  may  be  defined  by 


(199b) 


k  = 


(200) 


and  the  four  roots  of  the  characteristic  equation  are  k,  k,  — k, 

and  -k,  where  a  bar  indicates  the  taking  of  conjugates.  The  solution 
of  equation  (196)  may  be  written 

Ox  =  Cie-kx  +  T^e-*™  +  C2ekX  + 


(201) 
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Where  equation  (199)  does  not  hold,  which  is  the  case  for  very  small 
values  of  Gc/Ef  only,  all  four  roots  of  equation  (198a)  are  real  and  of 

the  form 


kl  = 

\|m l2  +  ' 

k2  = 

-ki 

k3  = 

\lm l2  "  ' 

m^  —  bm^ 

k4  = 

~k3 

(202) 


and  the  solution  of  equation  (196)  can  he  taken  in  the  form 


kix  — kix 

Qyr  —  Aje  +  Age 


k?x 

+  A^e  ^ 


— kox 

+  A4©  J 


(203) 


Before  applying  either  solution  to  a  specific  problem,  there  are 
noted  the  following  relations  which  follow  from  equation  (200) 


kk  =  | k| ^  =  2m2^ 

k  +  k  =  \J~2  \j  mx2  +  2m22 


(204) 


Seml-in-r-tnlte  shell  acted  upon  hy  edge  bending  moment  and  shear 
force.—  There  are  the  following  boundary  conditions. 


Mx(0)  =  &  +  *■)?*--  Pac*  (0)  =  Mq 
1  +  X  -  V2 


(205) 


Qx(0)  =  Qo 


while  for  x  =  «>  these  same  quantities  vanish. ^5 

15For  the  same  problem  without  the  effect  of  transverse  shear  and 
normal  stress^  see  reference  9* 
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Of  particular  interest  in  this  solution  are  the  values  of 
deflection  w(0)  and  change  of  slope  Px(0)  at  the  section  vhere 

the  loads  M, q  and  Qq  are  applied. 

Taking  first  the  case  Ef/Gc  <  2a/t  for  vhich  equation  (20l) 
applies j  it  is  seen  that  the  conditions  at  infinity  require  that 


C2  =  Cr}  as  0 


(206) 


so  that 


— kx  _  — kx 

Qx  =  ^1®  +  Cj^e 


(207) 


The  values  of  may  "be  obtained  "by  integration  from  equation  (l80) 

in  the  form 


(1  +  X)D*  fi  _  Ci  tjc  Ei  ^ 

,  o  px  ~  ,0  e  +  pe 

1  +  \  -  V2  kd  ^ 


(208) 


where  two  constants  of  integration  have  been  discarded  to  satisfy  again 
the  conditions  at  infinity. 

With  equations  (207)  and  (208)  there  is  obtained  from  the  boundary 
conditions  (equations  (205))  that 


Ci  +  Cq  =  Qo 

(Ci/k)  +  (c2/kj 


(209) 


•'■^Without  transverse  shear  and  normal  stress  deformation  these 
relations  are 


v(0)  -  \fw 


V2 


Ni  D*C* 


4(l-  v2)c* 


in  agreement  with  equations  (236)  of  reference  9 ,  where  the  homogeneous 
shell  is  considered. 
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This  determines  C-j_  and  in  the  form 

C  kQ0  +  1;12mo 

^  k  —  k 

-  kQo  +  |k|2Mo 

ci  =  — = — : — 


(210) 


Equation  (210)  is  introduced  into  equation  (208)  and  there  is 
obtained  as  the  first  of  two  "influence  coefficient"  formulas 


f1  + . X-)D-*P  3X(0)  =  - 

1  +  x  -  V2 


-Jpfoo  +  (k  +  k)Mo] 


(211) 


The  second  of  these  formulas  follows  from  equations  (l84),  (207) } 
and  (210)  in  the  form 

_^L_vM  =  _[jk|2Mo+  (k  +  k)QoJ  (212) 


Equations  (21l)  and  (212)  may  he  written  in  more  explicit  form^ 
using  equations  (204)  and  (197) •  The  results  are 


MCA  TN  No.  1832 


71 


and 


Neglecting  the  generally  small  effect  of  finite  Ec  in 

equations  (213)  and  (214),  that  is,  putting  X  =  0  in  these  equations, 
there  may  he  written  instead 


Equations  (215)  and  (216)  contain  the  noteworthy  fact  that  the 
correction  factors  for  the  effect  of  transverse  shear  are  independent 
of  the  ratio  t/h  of  face— layer  thickness  to  core  thickness.  The 
complete  formulas  of  course  must  and  do  contain  the  influence  of  the 
core  thickness  h. 

It  is  further  noted  that,  while  equations  (211)  to  (2l6)  have  been 
derived  for  the  case  that  m^  <  4mrA,  for  which  the  complex  solution 

holds,  they  are  also  valid,  as  is  readily  shown,  when  >  m^. 

Comparing ‘equations  (213)  and  (214),  and  (215)  and  (216)  with  the 
equations  listed  in  footnote  1 6  it  is  seen  that:  (l)  The  effect  of 
transverse  shear  modifies  the  deflection  due  to  Qq  and  the  rotation 

due  to  Mo  hut  not  the  other  two  coefficients,  (2)  the  effect  of 
transverse  normal  stress  enters  all  four  coefficients  hut  only  in  a 
minor  way,  and  (3)  the  reciprocity  relation  that  the  deflection  due 
to  Mq  is  the  same  as  the  rotation  due  to  Qq  is  carried  over  from 

the  theory  without  the  extra  effects. 
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For  a  numerical  example  the  following  data  are  chosen:  t  =  0.1  inch, 
h  =  1  inch,  a  =  10  inches,  Ef/Ec  =100,  Ef/Gc  =  200,  V  =  1/3, 

and  X  =  ^  •L,1IqQQ-‘^  100  =  O.O55,  and,  from  equation  (197), 


_  JL  0.1  X  200~ 
10  \j 1.018.X  1.1 


0.426 


m  kf 
*2=\] 


100  X  1.21  X  1.018 


0.294 


Then,  according  to  equation  (204), 


|k|2  =  0.173 

k  +  k  =  \J~2  V  0.182  +  0.173  =  0.84 


while  without  transverse  shear  deformation  (m^  =  0)  the  value 

of  k  +  k  =  0.59.  According  to  equations  (211)  and  (212),  the  effect 
of  transverse  shear  in  this  case  is  to  increase  the  rotation  due  to 
the  edge  moment  in  the  ratio  0.84/0.59  =  1.42,  an  effect  of  42  percent. 
The  same  increase  is  found  for  the  deflection  due  to  the  edge  shear 
force.  Rotation  due  to  the  shear  force  and  deflection  due  to  the 
moments  are  practically  unchanged.  Likewise,  the  effect  of  transverse 
normal  stress  in  this  case  is  of  negligible  importance. 

As  a  further  numerical  example  there  is  chosen  t  =  0.05  inch, 
h  =  1  inch,  a  =  20  inches,  Ef/Ec  =  1000,  Ef/G-C  =  2000, 

and  X  =  —  —  ■— 1000  =  O.O65,  and,  from  equation  (197), 

2  400 


mi  =  ■  2000  =  0.483 

1  20  V  1.022  X  1.05 
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From  equation  (204)  then 


|k|2  =  0.095 


0.234  +  0.095  =  0.82 


while  without  transverse  shear  deformation  (m^  =  0)  the  value 

of  k  +  k  =  0.44.  Thus  the  effect  in  this  case  is  to  increase  edge 
rotation  due  to  edge  moment  and  edge  deflection  due  to  edge  shear 
force  in  the  ratio  0.82/0.44  =  I.87,  an  effect  of  87  percent. 

Infinite  circular  cylindrical  shell  acted  upon  by  transverse  line 
load.—  Calculation  is  restricted  to  the  determination  of  deflection  and 
"bending  moment  at  the  section  x  =  0  where  the  line  load  of  intensity  2Qo 
is  assumed  to  act.  The  result  of  the  foregoing  paragraph  may  "be  used  as 
follows.  Consider  the  infinite  shell  cut  in  two  parts  at  the  section  x  =  0 
and  assume  a  tending  moment  Mq  of  such  magnitude  that  the  slope  Px(0) 

is  zero.  According  to  equation  (21l),  this  giveB 


Mq  =  - 


Qo 


Qo 


k  +  k  V 2mx2  +  4m22 


(217) 


and  therewith 


SSsdQl  ,  +  l*|2  0„  ,  _  a»i2  +  a»22  ,JA 

(1  *  l/3).2  .  k  +  k  ^2  -  hm22 


(218) 


Equations  (217)  and  (2l8)  become,  with  equations  (204)  and  (197) y 


iO 


c* 


Qo 

1  +  X  -  v2  c*/a2 


1 


(219) 


7k 


HACA  TN  Wo.  1832 


w(0)  =  - 


(*  *  iY  \J 


c* 


/  X  +  A;  —  V 

\/  a2D*  (1  +  X)(l  +  X/3) 

(h  +  t)G0  1  +  X/3 

,r 

C*  1  +  X  -  V2 

,  C*/a2  1 

r\ 

a2D*  (1  +  X)(l  +  X/3) 

(h  +  t)Gc  1  +  V3 

To  give  these  formulas  a  somewhat  less  unwieldy  appearance,  the 
effect  of  finite  Ec,  that  is,  X  =  0,  may  again  he  neglected,  as  is 
permissible  in  most  cases;  and  there  may  he  written 


Mq  = 


w(0) 


\J  1  -  V2  a2 

^  t\jjh  +  t)a 


1  + 


1 

\/l  -  v2 


t  ^f 
a  Gc 


t  ff 
a  Gc 


Qp 

Ejp 


Some  numerical  examples  are  given  as  follows. 

Taking  t  =  0.1  inch,  a  =  10  inches,  Ef/Gc  =  200,  and  V  =  l/3, 
transverse  shear  deformation  reduces  Mq  to  1 / \/  2 . 05  times  the 
value  which  holds  when  Gc  =  00,  that  is,  there  is  about  a  30  percent 
reduction  in  Mq.  At  the  same  time  the  deflection  under  the  line 

load  is  3.05/\]2.05  =  2.lh  times  what  it  is  when  Gc  =  00;  that  is, 
there  is  an  increase  of  about  115  percent  in  w(0). 

Taking  t  =  0.05  inch,  a  =  20  inches,  and  Ef/Gc  =  200,  Mq  is 
decreased  by  a  factor  \fk/[ 5  =  0.89,  while  w(0)  is  increased  by  a 
factor  1.5/ v/l.25  =  1.3^-» 

Taking  t  =  0.05  inch,  a  =  20  inches,  and  Ef/Gc  =  2000,  Mq  is 

decreased  by  a  factor  l/sJJjfe.  =  O.526,  while  w(0)  is  increased  by 
a  factor  6. =  3.29. 


(220) 


(221)* 


(222)* 
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Equation  (220)  for  w(0)  may  be  compared  with  equation  (ll6)  for 
the  circular  plate  of  radius  a.  This  comparison  shows  that,  while  for 
the  plate  both  the  ratios  t/a  and  (h  +  t)/a  enter  into  the  correction 
factor,  the  correction  factor  for  the  cylindrical  shell  contains  the 
ratio  t/a  only;  that  is,  the  corrections  (but  not  the  results)  are 
independent  of  the  ratio  of  face-layer  thickness  to  core  thickness  in 
this  case  of  a  cylindrical  shell. 


Spherical  Shells 

In  conformity  with  customary  usage,  the  following  notation  is 
introduced : 


61  • 

a0 

£2  = 

a0 

al  = 

1 

02  = 

sin  0 

Ei  = 

E2  =  a 

Nn 

ll 

n22  ■ 

=  N 0 

N^  , 

■  W21  =  N00 

%  = 

00 

02  = 

% 

Mil  > 

=  M0 

%>  ■ 

=  «e 

%2  = 

-  ttye 

Pi  = 

P0 

P2  = 

P0 

ml  = 

y 

D32  = 

“0 

U1  = 

U 

u2  = 

V 

Pi  = 

P0 

P2  = 

Pe 

Attention  is  here  restricted  to  problems  with  rotational 
and  the  following  relations  are  used: 


(223) 


symmetry 


S(  )/be  =  0 

N0e  =  %  =  u<Ae  =  vq  =  =  v  =  p0  =  0 


(224) 


The  differential  equations  of  equilibrium  (25)  to  (28)  (31) 

and  (34)  become,  setting  *  3 


d(  )M  =  a(  )/d0  =  (  )* 
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There  is  first  given  a  simple  special  solution  of  this  system  of 
equations  and  then  a  generalization  is  obtained  of  the  two  simultaneous 
equations  for  and  (3^  which  are  fundamental  in  the  theory  of 

homogeneous  isotropic  Bhells. 

Uniform  stress  distribution  in  a  spherical  shell.—  Set  in 
equations  (225)  to  (23*0  P0  =  m^  =  0  and  assume  that  N^  Ng,  Qty, 

»nii  Mg  are  independent  of  $ .  From  equation  (225)  it  follows  that: 


(235) 


From  equation  (226)  it  follows  then  that 

Nq  =  |  aq  (236) 

and  from  equation  (227)  it  follows  that 


M^  =  Mg  =  Mq 


(237) 


Equation  (228)  gives 


=  s  -  2Mo/(h  +  t)a 


(238) 


In  equations  (230)  and  (23l)  set  u  =  0  for  reasons  of  symmetry  and 
obtain 


(l  +  §X-v)N0. 


c* 


W  + 


t(h 


6a 


+  t)  Ef 
E„ 


(239) 


I 
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or,  with  N0  from  equation  (236)  and  \  from  equation  (229), 


(C*/a)w  =  |(l  -  V)aq 


(240) 


Equation  (232)  is  identically  satisfied  when  (30  =  0.  Equations  (233) 
and  (234),  in  conjunction  with  equation  (237),  give 


(1  +  24  -  V)Mq  =  (D*/a)(s/Ec) 


1  t(h  +  t)2  Ef 

2  a  ec 


or 


M 0  = 


(h  +  t)a4 
1  +  24  -  v 


s 


(241) 


Then,  from  equation  (238), 


°tm  = 


...  (1  - 

1  +  24  -  V 


(242)  , 


Equation  (242)  may  he  compared  with  equation  (136)  for  the  circular  ring. 

According  to  equation  (46),  there  are  obtained  from  equations  (236) 
and  (24l)  the  following  expressions  for  the  stress  resultants  in  the 
outer  ("upper”)  and  inner  ("lower")  face  layers: 


(243) 


Comparison  of  these  results  with  the  corresponding  results  for  the 
circular  ring  (equations  (135))  shows  that  for  given  values  of  q 
and  s  there  is  a  greater  difference  between  Nu  and  N^  in  the 
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spherical  shell  than  there  is  in  the  circular  ring,  the  reason  "being 
the  relatively  larger  influence  of  the  s-term  in  equation  (243). 

For  a  specific  example,  it  is  again  assumed  that  the  radial  load  is 
applied  entirely  to  the  inner  face  so  that  qu  =  0  and,  according  to 

equations  (20)' and  (22), 


Substitution  of  equation  (244)  in  equation  (243)  gives 


h  +  t 
2a 


H  ■  (> 


h  +  t\2 
2a  ) 


a  <Li 

T 


1  -  V 

1  +  2X  -  V 


h  +  t\a<U  1  +  4l  -  v 
2a  J  k  1  +  2X  -  v 


(245) 


As  a  numerical  example,  taking  X  =  O.O595,  as  in  the  example  given  in 
the  section  entitled  "Closed  circular  ring  acted  upon  by  uniform  radial  load 
and  V  =  1/3,  it  is  found  that  the  factor  in  Nj  which  contains  the  effect 
of  the  core  flexibility  is  (l  +  0.36)/(l  +  0.18)  =  I.15.  Thus,  where 
for  the  circular  ring  there  was  a  6rpercent  stress  increase,  there  now 
is  a  15-percent  stress  increase. 

Reduction  of  axisymmetrical  problem  to  two  simultaneous  equations 
for  Q0  and  P0.—  The  fundamental  results  of  reference  10  for 

homogeneous  shells  may  be  readily  extended  to  sandwich  shells  as 
follows : 

Equations  (225)  and  (226)  are  used  to  express  N^  »nri  N0  in 
terms  of  Q 0. 

w0  =  cot  0Q0  +  Fx(0)  (246) 


N0  =  Q0*  +  F  2(0) 


(247) 
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In  equations  (246)  and  (24?)  the  functions  Fj  and  F2  are  given  by 


Fi  = 


sir 5r0 


J*  (q  cos  0  —  p0  sin  0^sin  0  d0 


(248) 


*2  " 


cos  0 


Jjsin  0  Fxy  +  ap^  sin  0j 


(249) 


Next  the  displacement  components  u  and  w  are  expressed  in  terms 
of  Q a,  "by  means  of  equations  (230),  (23I) ,  (246),  and  (247)’. 


Subtraction  of  equation  (231)  from  equation  (230)  gives 


§V  -  cot  0  u)  =  (l  +  v)(N0  -  N@) 


=  (1  +  v)[-(Q0f  -  cot  00^)  +  Fi  -  Fg|  (250) 


Equation  (112)  is  integrated  to 


(C*/a)u  =  -(1  +  v)(Q0  +  F3) 


(251) 


where  Fo  is  given  "by 


P  Fx(0)  -  F 2(0) 

Fo  =  —sin  0  - - -  d0 

3  r J  sin  0 


(252) 


Equations  (25l)  and  (252)  are  introduced  into  equation  (23l)  and  the 
following  expression  is  obtained  for  w: 


(C*/a)w  =  (1  +  X/3)^cot  0Q0  +  Q0f^  +  F^ 


(253) 
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where  F^  is  given  "by 


Fi.  =  - 


■(h3S-^  ^  +  (1  +  v)cot  0  f3  +  (l  +  §)F2  -  (v  -  |)f- 


(254) 


Equations  (251)  and  (253)  are  introduced  into  equation  (232)  for 
and  the  first  of  the  two  simultaneous  equations  for  and  00  is 

obtained  in  the  form 


(h'7t)0c  -  +  -  f)  (cot  <¥ +  y)’ +  * 

+  (1  +  v)  (00  +  f3) 

which  may  h©  rearranged  to  read-*- 7 


Q0» »  +  cot  0  Q0*  -  [< 


cot20  -  1-  X/3  +  SI _ 1_ 

1+1/3  Gc  (h  +  t)(l 


(1  +  1/3)] 


1  +  l/3 


P0  =  F5(0) 


(255) 


The  function  F^  is  given  hy 


w  C1  +  V>F3  +  F1 

5  ■  —  ;  x/3 — 


(256) 


Introducing  the  operator 


L  =  (  ) ’  +  cot  0(  ) ’  -  cot20(  ) 

ITWhen  1=0  and  Gc  =  00  and  when  no  external  loads  are  present^ 
this  equation  checks  with  the  first  of  equations  (g)  on  page  t-69  of 
reference  9. 
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Again,  using  the  operator  L,  this  may  he  written  in  the  form 

(l  -  p[l  -  V2  +  2X(i  +  V)]Q0  =  Fg(0)  (262)* 

The  function  Fg  is  given  hy 

Fg  =  -(1  +  V)f£  -  |§[l  —  V2  +  2X(1  +  v)]ny  (263) 

Equation  (262)  may  he  compared  with  the  second  of  equations  (g) 
on  page  ^69  of  reference  9. 

Analysis  of  edge  effect  for  spherical  shell.—  The  special  case  of 
no  distributed  surface  load  and  no  concentrated  load  at  the  apex  of  the 
shell  is  obtained  hy  setting 


f5  =  Fg  =  0 

Following  again  a  known  procedure  from  the  theory  without 
transverse  stress  deformation,  there  may  he  set 
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with  corresponding  formulas  for  Pgrf  and  P /  .  Introduction  of 

equation  (264)  into  equations  (257)  and  (262)  gives 

%”  +  (m  +  §  ~  |  cot2^)%  +  3-7773  h  =  0  (266) 


Pl. »  _  -1  +  2  cot20jP!  -  j|[l  -  V2  +  24(1  +  v)]  Qq  =  0  (267) 


Attention  is  restricted  to  the  cases  for  which  =  Fg  =  0  in 
equations  (266)  and  (267)  and  the  problem  is  considered  of  the  shell 
subject  to  edge  loadB  (^0)oJ  an<^  (^0)o  a  sec^lon  9  ~ 

Assuming  that  cot  0O  is  no1:;  larB©  compared  with  unity1®  and  that  the 
effect  of  the  edge  loads  is  restricted  to  a  narrow  edge  zone  so 

that  |  Qx|  «  |Qi»’|,  I  Pd|  «  |p1,,|,  equations  (266)  and  (267)  may  be 

simplified  to 


%"  -  +  T^rz  h  -  0 

1  +  4/3 


(268)* 


+  24(1  + 


0 


(269)* 


Equations  (268)  and  (269)  show  that  the  influence  of  finite  Ec(4  /  0) 
in  the  edge-effect  problem  consists,  except  in  extreme  circumstances, 
in  minor  modifications  of  the  results  for  Ec  =  00.  The  quantity 

which  represents  the  influence  of  finite  Gc  and  which  is 


~  2t  Ef 
P1  '  h  +  t  ttc 


(270) 


xna y,  however ,  in  practical  cases  he  large  compared  with  unity  and  not  of: 
negligible  influence  on  the  results. 

-^Afhen  cot  »  1  the  shell  is  termed  a  "  shallow”  shell  which  is 
not  considered  in  what  follows. 
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,  ,E|u®tions  (2^8)  and  (269)  may  "be  compared  with  equations  (l8o) 
and  (182)  for  the  cylindrical  shell.  This  comparison  shows  that  the 
influence  of  finite  Gc  in  the  edge-effect  problem  is  of  the  same 
-nature  for  the  spherical  and  cylindrical  shells.  Thus,  results  of  the 
same  quantitative  nature  will  he  obtainable  as  in  the  section  on 
cylindrical  shells  under  the  headings  entitled  "Finite  circular  cylindrical 
shell  acted  upon  by  edge  moments  and  forces"  and  "Semi-infinite  shell  acted 
upon  "by  ©dg©  *b ending  moment  and  shear  force.” 

This  work  is  not  herein  carried  further  to  specific  applications. 

SUvh  aPlications  “T  136  worked  out  with  hardly  any 
more  difficulty  than  when  the  effect  of  the  core  deformability  is  not 
taken  into  account. 


CONCLUDING  REMARKS 


A  system  of  basic  equations  has  been  derived  for  the  analysis  of 
small-deflection  problems  of  sandwich-type  thin  shells.  This  system  of 
equations  reduces  to  Love*s  theory  of  thin  shells  when  the  transverse 
shear  and  normal  stress  deformability  of  the  core  of  the  sandwich  is  of 
negligible  importance.  The  system  of  basic  equations  has  been  applied 
to  a  number  of  specific  problems  from  the  theory  of  plates,  circular 
rings  circular  cylindrical  shells,  and  spherical  shells,  and  it  has 
been  found  that  the  effects  of  both  transverse  shear  and  transverse 
n°™f\pr08a  Reformation  may  be  of  such  magnitude  that  an  analysis 
which  disregards  them  gives  values  for  deflections  and  stresses  which 
are  appreciably  in  error.  . 

Numerical  calculations  have  been  in  the  nature  of  sample  calcu- 
lations,  illustrating  both  the  use  of  the  equations  and  the  possible 
effects  of  using  them.  Examples  have  been  choBen  from  the  point  of 
view  of  relative  simplicity  as  well  as  with  the  thought  to  illustrate 
most  clearly  the  consequences  of  the  extra  deformations  which  have 
been  taken  into  account.  It  is  unavoidable  that,  in  so  doing,  some  of 
the  examples  may  be  of  little  interest  for  aircraft  structural  analysis 
and  that  some  problems  may  not  have  been  analyzed  which  would  have  well 
fitted^within  the  contents  of  this  report  and  which  at  the  same  time 
would  have  been  of  considerable  practical  importance. 

The  general  analysis  has  been  restricted  by  the  following  two 
order-of-magnitude  relations:  (l)  t/h  «  1  and  (2)  tEf/hEc  »  i} 

where  t  is  the  face-layer  thickness,  h  is  the  core-layer  thickness, 

■&f  is  the  elastic  modulus  of  the  isotropic  face-layer  material 
and  Ec  is  the  elastic  modulus  in  the  transverse  direction  of  the 
core-layer  material  Therewith  it  is  felt  that  very  likely  nearly 
all  situations  have  been  covered  in  which  the  effect  of  transverse 
core  flexibility  is  of  significant  practical  importance.  It  is 
evident,  however,  that  if  desired  the  theory  could  be  extended  so  as 
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Jfc 


llfle  caSe3  where  one  or  both  of  these  two  order-of-Hnagnitude 
™JtiS  are  not  satisfied .  The  min  imitation  of  the  P^sent 
^rmthe  omiBsion  of  all  finit^eflection  and  instability 

effects. 


Massachusetts  Institute  of  Technology 
Cambridge,  Mass.,  May  26,  19^7 
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